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(— I , Abstract 

We give a detailed study of dynamical properties of the Zhang model, 

S including evaluation of topological entropy and estimates for the Lya- 
punov exponents and the dimension of the attractor. In the thermody- 
namic limit the entropy goes to zero and the Lyapunov spectrum collapses! 

m ■ 
> ' 

O ! Introduction 

G\ ■ 

' In 1987 the concept of Self-Organized Criticality (S0C) was introduced by Bak, 

Tang and Wiesenfeld |BTWj . The attempt was to give an explanation of the 
omnipresence of fractal structures and power-law statistics in nature, and the 
claim was that certain physical systems can self-organize into stationary states, 
reminiscent of equilibrium system at the critical point, in the sense that one has 
scale invariance and long range correlations in space and time. 

S0C is proposed as an explanation for variety of phenomena in nature, such 
as earthquakes, forest fires, stock markets and biological evolution PJ. However, 
most work has been devoted to the study of idealized " sandpile-like" computer 
•*h . models, such as the sandpile model |BTW] . the abelian sandpile |DR| and the 

Zhang model [7] that, one believes, exhibit S0C in the thermodynamic limit. 
Despite this effort, a satisfactory understanding of the model is not yet achieved. 
Through numerical investigation it was observed that in the thermodynamic 
limit, observables have power-law distributions. More precisely, the probability 
distribution of an observable s has the form P(s) ~ 1/ s Ts in the thermodynamic 
limit. There is no widely agreed upon method for computing the SOC-exponents 
t s numerically and, due to the incomplete understanding of the dynamics of the 
models and lack of a formal treatment of the thermodynamic limit, it is difficult 
to properly explain the observed behavior. Hence it is not clear what the S0C- 
cxponcnts really tell us about the dynamics of the S0C models. 

1 Keywords: sand-pile models, avalanche dynamics, skew- product systems, Lyapunov ex- 
ponents, entropy, HausdorfT dimension, thermodynamic limit. 
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0.1. Discussion of the Zhang model 

In a series of papers by Cessac, Blanchard and Kriiger |BCK| it was proposed 
that deeper understanding of SOC models can be achieved by studying the 
models in the framework of dynamical system theory. They showed how a 
particular model, the Zhang model, could be formulated as a dynamical system 
of skew-product type with singularities, where the randomness of the external 
driving is described by a Bernoulli shift, and the threshold relaxation dynamics 
is given by piecewise afhne maps. 

In this paper we present a detailed study of the dynamical system defined 
in |BCK| . We prove several basic properties, some of which are already stated 
in |BCK| . before discussing fundamental dynamical properties. Depending on 
the parameters of the model, we can observe fundamentally different types of 
behavior. 

For low values of the threshold energy (critical energy), the dynamics can 
be relatively simple, since the singularities only effect the dynamics in a finite 
number of time-steps. In such situations we say that singularities are remov- 
able, and we show that the system permits symbolic coding. We give examples 
of how symbolic coding provides a complete description of the dynamics as a 
topological Markov chain. Hence the dynamics is chaotic, but the essential dy- 
namical invariants are all inherited from the Bernoulli shift factor. Moreover 
we can identify the physical invariant measure, and hence understanding of the 
statistical properties is reduced to the theory of Markov chains. 

As we increase the critical energy the role of the singularities becomes essen- 
tial. Techniques based on codings are no longer applicable and a very interesting 
dynamics emerges. The dimensional characteristics of the attractor are also sen- 
sitive to the parameters of the system, as we show by generalizing the Moran 
formula for the iterated function system (IFS). In addition, we observe the situ- 
ation, when the dimension of the IFS-attractor increases to the maximum, while 
the support of the SRB-measure remains fractal. 

To measure the complexity of the dynamics we study entropy and Lyapunov 
exponents. We show that the system is hyperbolic, with one positive exponent 
originating in the Bernoulli shift. However, due to the presence of singularities 
the Ruelle inequality and the Pesin formula are not directly applicable. We 
show that the metric entropy of any SRB-measure equals the topological entropy 
almost surely, and we evaluate the latter generalizing the technique developed 
by Buzzi |B1I IB2| . The result is that the Pesin formula and the variational 
principle hold a posteriori. 

To give a satisfactory physical interpretation of the dynamics we rescale time 
to prevent infinitely slow driving of the system. We prove that for this physical 
system, the Lyapunov spectrum collapses completely and that the entropy goes 
to zero in the thermodynamic limit. This implies that the expanding (chaotic) 
properties are lost, so that we may expect power-laws statistics and long range 
correlation effects. 

The statistical properties we obtain hold for any SRB-measure, because the 
most input comes from the Bernoulli shifts. The existence of SRB-measures is 
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in fact still an open problem. From the general theory of dynamical systems 
with singularities |KSMPlllST| . we can give conditions that are sufficient for the 
existence of SRB-measures, but it is not known if these conditions hold for the 
majority of parameters. We expect this to be true (it was also conjectured in 
jBCKp and derive some statistical corollaries. 

Apart from the physical importance of the Zhang model, it is interesting 
from a mathematical point of view. It can be described as a piecewise affine 
hyperbolic map of the form 

F : E+ x M -> E+ x M, {{t ht 2 ...),x)~ ((*ii 2 t 3 ■ • • ), /*„ (a:)), 

where E^- is the set of right infinite sequences from a finite alphabet and {/,} 
a collection of piece-wise affine non-expanding maps of M to itself. Previously 
piecewise affine expanding maps and piecewise isometries have been studied, 
but the contracting property of the relaxation dynamics gives rise to some dif- 
ficulties. Therefore several methods are developed in this paper, which hold far 
beyond the framework of the Zhang model. 

0.2. Structure of the paper 

In Section ^ w e describe the model and derive bounds on the size and duration 
of avalanches. This enables us to use the Poincare return to reformulate the 
systems in a skew-product form. Then we study the contraction property to 
conclude hyperbolicity of the model (Theorem and describe, when degen- 
erations occur (Theorem [SJ the original Zhang setting e = is not the only 
possibility). In Section [3 we introduce the concept of removability of singulari- 
ties, which appears in the coding approach for the study of the model. 

Section |21 is devoted to the study of measure entropy and Lyapunov spec- 
trum. We prove in Theorem 1 1 31 that the entropy of an SRB-measure is always 
maximal. Section 0] concerns the topological entropy. We evaluate it for the 
most parameter values (Theorems El an d I19|l . We also discuss nearly-Zhang 
models and show that the dynamical quantities do not change. This is natural 
from the physical perspective, because SOC should not be obtained through a 
fine tuning of parameters. 

Section [3] briefly describes the dimension issues of the model (Theorem EU1 
gives the asymptotic values), which enters into all inter-relations involving en- 
tropy and characteristic exponents. We demonstrate how the fractality occurs in 
IFS-context, noting the difference due to singularities and overlaps. In Section 
El we illustrate the most important effects in the model by examples. 

In Section we discuss the thermodynamic limit and attempt to explain 
appearance of the power-law statistics via a reparametrization. Conclusion con- 
tains the physical implications of the current investigation. 

In Appendices[^and[B]we provide bounds for the entropy and the dimension, 
which are new in the presence of singularities, overlaps and degenerations (this 
was designed for an application to the Zhang model) . The results are of interest 
in its own and can be read independently. 
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1. Basic properties of the Zhang model 

In the Zhang model each site on the lattice is associated with a non-negative real 
number, which we call the energy of the site. The collection of energies is called 
an energy configuration, and can be represented as a point in ./V-dimensional 
space, where N is the number of sites in the lattice. If a configuration is unstable, 
the overcritical sites will lose some of their energy to their nearest neighbors, 
resulting in a new energy configuration. This transformation on is denoted 
by /. If a configuration is stable, a site is chosen at random and an energy 
quantum S — 1 is added to this site. In |BCK) it was shown how the relaxation 
and random excitation can be formulated as a map of skew-product type on an 
extended phase-space. This extended phase-space has the configuration space as 
one factor, and the set of all possible sequences of excitations as the other factor. 
In |BCK| it was also shown how one can reformulate the dynamical system by 
considering the return maps to the set of stable configurations. This gives 
a simplification, in the sense that each avalanche is associated with an affinc 
transformation. The set of stable configurations is partitioned into domains, 
where each domain corresponds to an avalanche. 

1.1. Relaxation 

Take d,LeN and let A C 1 d be the cube [1, L] d of cardinality N := L d = |A|. 
Let c/> '■ A — * A' := {1, . . . , N} be a bijection. We define a metric d\ on A by 



and let d^ :— (f>^d^. In the following we omit primes when it is clear from the 
context that we are considering the metric space (A',d\>). Elements of A will 
be called sites. We say that sites i and j are nearest neighbors if d\(i,j) = 1. 
The boundary dA is defined as those sites i £ A that have less than 2d nearest 
neighbors. 

Fix parameters E c > and e 6 [0, 1) and define / : R> — > R> by 



<iA(k,l)= 2_j \ k ri~ln\ 



Kn<d 



Xi - 6(xi - E c )(l 



1 - e 



E 



6(x 



E c )x, 



2d 



d.A(i,j) 



1 



where 




Let ||x||i = Yli=i \ x i\ b e the 1-norm on R 
Proposition 1. For all x £ R> we have 



1 + e 
2 



N|i < ||/(as)||i < Hi, 
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and ||/(.t)||i = ||x||i if and only if Xi < E c for all i € OA. If there is i G <9A 
such that Xi > E c then 

ll/(»)Hi< INIi-^^c 

Proof. Let {xi k }™ =1 be the entries of the vector x that are greater than E c . 
Let rii k be the number of nearest neighbors of Xi k . Then 

m _^ m 

ieA ieA k=l k=l 



= ^-(l-e^l-^K 



2d 

ieA fe=i 

The statement follows from the fact that we always have d < n ik < 2d, and 
n ik = 2d if and only if Xi k £ dA. □ 

We say that a site i S A of the configuration x is relaxed if Xi < E c , and 
excited if xi > E c . A configuration x is called stable if all sites are relaxed. The 
set of stable configurations is M := [0, E c ] . For each configuration x we define 
m{x) = mm{n > | f n (x) € M}. 

Proposition 2. For all x e M.> we have: 

, s 2dA^ ||ic||l / 2d xdiam(A)/2 

We need the following lemma ([•] denotes integer part): 

Lemma 3. For x £ M> and n G N let ai(n,x) be the cardinality of the set 
{I < n | (fx)i > E c }. Let 7 = [2d/(l - e)] + 1. Ifd A (iJ) = 1, then a 3 (n,x) > 
[ai(n, x)/i\. 

Proof. There is a finite increasing sequence {mk}, such that (f mk (x))i > E c . 
We claim that on each interval {mk 1 mk+ 1 \ there is a number m such that 
(f m (x))j > E c . In fact, in the opposite case 

( / l+m !!+ ,-x (a . ) ).> 7 i__f Sc>jE ; c . 

Since [0, ai(n, x)] contains [3 — [cti(n,x)/"f] disjoined such intervals, we get 
ctj(n,x) > [3. Thus aj(n,x) > [ai(n,x)/j]. □ 

Proof of Proposition [2J By applying inductively Lcmma[3]we get: 

"aj(n, x) ' 



aj(n, x) > 
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In fact, if j = jo, ji, ■ ■ -jk = i is a path with d A (js Js+i) = 1 and 

-ai(n,x)- 



t. 



then aj k [n, x) > , aj k _ 1 (n, a;) > t-f k 1 , . . . , ctj (n, x) > t. By Proposition^ 

/ w 2d ||x||i 

for j 6 9A, so 

for all i e A and all neN. Suppose f m (x) £ M for all m < T. If T > iVa(x), 
then there must be a site i £ A that is greater than E c for more than a(x) 
different times. This is impossible so m{x) < Na(x). □ 

1.2. Random excitations 

Define E^- = A N to be the set of right-infinite A-sequences and let : E^- — ► 
E+ be the left shift. We define a map / : E+ X R^ -> E+ x M^ by 



/(M) = 



(o^t, a; + e to ) if x e M 
(t,/(ar)) i£x#M' 



where ei, . . . , is the standard basis in M. N . We denote points in E^- x R> 
by x — (t,x), and we define tt u and 7r s to be the projections to E^- and R> 
respectively. 

Proposition 4. For all £ € E^ x R> holds: 

mm{m > | Vi S A 3m' < m : (tt s o / m ' (£)), > £? c } < n(.E C) e, A) , 

where 

T 2d 1 x diam(A) 



n(£ OJ e, A) = N(NE C + 2) ( [— ] + I ) 



Proof. In N[E C ] + I time-steps, there must be an overcritical site. Since in 
the relaxation process there is always an overcritical site, then during arbi- 
trary subsequent N[E C ] + 2 time-steps an exited site can be found. Hence after 
N£(N[E C ] + 2) time-steps either all sites have been overcritical or there is a site 
that has been overcritical at least £ times. However it follows from the proof of 
Proposition |2 that if one site is overcritical 

2d 



diam(A) 



f 



.1 -e- 

times, then all sites have been overcritical at least once. □ 
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For x G M> and i G A we define r(i, x) := min{n € N | f n (x + ej) G M}. 
Proposition |2 assures us that this number is finite and 

/ 1 \ / 2d \diam(A)/2+l 

maxmaxT(!,i) < T m (E c ,e,A) = N 2 [l + — — ( + 1 

x£M ieA V N E C J VI — £ / 

Thus we observe that neither n(E c , e, A) nor r m {E c , e, A) are uniformly bounded 
in E c , but there is the following alternative: 

There exists a constant Co, not depending on the energy E c , such that either 
n(E c ,e,A) < Co or T m (E c ,e,A) < Co- 

In fact, we can set C = 3iV 2 (^ + !) diam ( A ) +1 . Thus we get that either 
relaxation happen sufficiently fast or all the sites keep being excited sufficiently 
often (uniformly in E c ). 

But there does not exist such a bound uniform in e or N. 

1.3. Return maps 

Let x = (t, x) G x M. For n = 1, . . . , T(t Q ,x) define C n (x) = {i G A | (tt s o 
f n x)i > i? c }, and A(x) = (C\(x), . . . , C T ( totX )(x)) . We call A(x) the avalanche 
of the point x. Let M := x M and define an equivalence relation ~ on M 

by 

x ~ y «^ = . 

This gives a partition of M. From the definition it is clear that A{x) depends 
on to and x only. Hence partition elements are of the form [i] x My , where 

Vi G A : |J My = A/ 

j 

and [i] = {t G | £o = is the cylinder of the symbol i. We see that for 
each i G A, the domains Mii,M i2 , . . . are separated by segments of at most 
iV! Tm = cxp (T m (E c , e, A) log Nl) hyperplanes. Hence we have a finite number of 
domains Mn, . . . , Mi qt for each i G A. By definition there is a unique avalanche 
for each partition element [i] x My. We denote this avalanche by Ay. Its 
duration is ry := r(i, x), for a; G My, and define its size to be sy = 53n=i |Cn|- 
We define the piecewise continuous map F : M — ► M by 

(t,x) ^ (cr^t,F to a;) 

where := f T ^ x \x + e,). We define Fy := F^My ■ 

Remark 1. From a mathematical point of view the formulation (M, F) is a 
simplification compared to (S^ x R> ,/). However, the duration of avalanches 
are suppressed so that all avalanches have the same duration. This is not satis- 
factory from a physical point of view, and hence we call (M, F) the mathematical 
model and (S^ x M> ' /) ^ e physical model. We will later make a rescaling of 
time in the physical model, so that the driving does not become infinitely slow 
in the thermodynamic limit. 
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For each x G M> we define a matrix Q(x) by 

Q kl (x) = l^ d{Xl ~ Ec) iid ^k,l) = l, 
1 otherwise, 

and a diagonal matrix J(x) by J k i{x) = (1 — (1 — e)9(xi — E c ))8 k i. Set 

S(x) = J(x) + (1 - e)Q(x) (1) 

and observe that /(x) = S(x)x. Let x(l) = x + and x(n) = f(x(n — 1)) for 
n G {2, . . . , r(t, a;)}. Then Fi(x) = Li(x + e;), where 

Li(a;) = S(a;(r(a;,t)))...5(a;(l)). 

If x,y G My, then r(i,x) = r(i,y) and the same components of x(n) and y(n) 
are grater than B c for each n = 1, . . . , r(t, x), so Lj(x) = Li(y). We define the 
linear map Lij := Li(x) for x G My. We get i^My ( x ) — Lij(x + e,). 

Definition 1. ^4 sequence {(i n ,j n ) \1 < n < 9} is said to be admissible if 
e 

n=l 

Theorem 5. For a/Z < 1. Moreover for every constant c G (0,1) 

i/iere is a number T G N swc/i i/iai /or every 9 > T and admissible sequence 
{{in,jn) 1 1 < Ti < 9} it holds: 

\\Lifjt ■ ■ ■Li 1 j 1 \\x < c. 

Proof. If A is an N x TV matrix we let C k (A) be its fc-th column. Observe that 
for any matrices A and B we have the following formula: 

WCkiAB)^ = 53||C,(i4)|| 1 B Ifc . (2) 
I 

By the construction: ||Cfc(S(x))||i < 1 for all k G A. Hence < 1. 

To prove the second statement we note that for e > the diagonal elements 
of the matrices S(x) are non-zero and > e. Therefore 

{S{x(m))S{x(m - l)))u > e • max{S , fei (x(m)), S k i(x(m - 1))}. 

Moreover, Ski{x(m)) > if x(m)i > and c£a(&, I) = 1. It follows that any ad- 
missible product Li e j e . . . Li 1 j 1 of length 9 > n(E c ,e,A) is positive. By Propo- 
sition n there must be at least one column such that the sum over this column 
is less than 1, for some factor £,: t j t and hence for the whole product. Therefore 
the sum over each column of any admissible product of length 2n(E c , e, A) must 
be less than 1. Let Co < 1 be the maximal norm of all admissible products of 
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length 2n{E c , e, A). For k > ko :— [logc/logco] + 1 we have c§ < c and hence 
T = 2kon(E c , e, A) is the required number. 

The above argument does not apply to the case e = 0, and a different proof 
must be given for this case (which actually works in general as well). Take 
x £ M and let x(t) £ M be the projection of its orbit to M. Denote S(x(t)) by 
St(x), and let 

S t (x) = S t {x) ■ ■ ■ S {x) . 
We make the following claims: 

1. There exists neN such that for all Z, m £ A and all x £ M there is t < n 
such that (St(x))i m 0- 

2. For all i > there exists n, £ N such that ||C m (S't(a;))|ji < 1 for all t > rii, 
x £ M and all sites m £ A with d\(m, dA) < i. 

The second claim for i — \ diam(A) implies the statement of the theorem. 

To see the first claim we fix x and let U C A 2 be the subset of the pairs 
(/, m) with (S t )im — for all sufficiently large t. By inductively applying @ we 
see that the columns for §t(x) are non-zero for all t > 0. So for all (3 £ A there 
is a £ A such that (a, (3) £ A 2 \ U. Given sites a and (3 we choose t such that 
St(x) a 7^ 0. Consider now column a of the matrix St+i(x) = St+i(x)St(x). If 
a is stable, i.e. x(t + l) a < E c , then (St+i(x)) aa > and (St+i{x)) a p ^ 0, so 
we just repeat the argument. But the site a can not be stable for more than 
n(E c ,Q,A) iterations. Hence we can with no loss of generality choose t such 
that x(t + \) a > E c . Then the column a of St+\{x) has non-zero elements in all 
position that correspond to neighbors of a. Hence we obtain that (a' , (3) £ A 2 \U 
for all a' with 4(tt',ct) — 1- Any two points can be connected by a path of 
neighbors, so U = 0, and the first claim follows. In fact, one can see that the 
bound n does not depend on a choice of x and satisfies: h < diam(A) -n{E c , 0, A). 

To prove the second claim let us note that if \\Ck{St{x))\\\ < 1, then 
||C fc (St+i(x))||i < 1 because by ©: \\C k {AB)\\ < max, \\Q(A)\\ ■ \\C k (B)\\. 

We will use induction on i starting from i = 0. Take k £ dA and t < n such 
that {S t (x)) kk + 0. If x(t+ l) k > E c , then ||C fc (5 t+1 (x))||i < 1 and 

||C fc (S t+1 (f))||i = \\Ci(S t+1 (x))US t (x)) lk < 1 , 
i 

and so we have the desired inequality. If x(t + l) k < E c , then (St+i(x))kk = 1 
and hence (St+i(x))kk ^ 0. Then we repeat the argument. Since no site can be 
stable for more than n(E c ,0,A) successive time-steps we obtain the claim for 
i = with tiq = n(E c , 0, A) + n. 

Consider now the case i > 0. For a site m £ A with d(m, dA) — i, we take 
I £ A with d\(l, m) — 1 and d\(l, dA) = i — 1. By the first claim we find some 
t < n such that (St{x))i m ^ 0, and by the induction hypothesis for t' > n,_i 
we have: 

||C,(S , t +t ,(f)-..5 t+ i(:c))||i < 1. 
Using (J5J) we obtain \\C m [St+t'(x )) ||i < 1. We can choose rii = ^i-i + n. □ 
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Lemma 6. Let E c > e/(l — e). Then for any x G M, n G N and i,j G C„(x) 
we /iat>e d\(i,j) 7^ 1. 

Proof. Take x G M and let £"„ be the maximal energy of a site in C n (x). 
Clearly E x < E c + 1 and 

E n +i < max { max{e£ n , £ c } + (1 - e)E n , E c + l} . 

From this we see by induction that 

E n < max { — ,E C + 1\ = — , 

so ei?„ < E c for all n G N, and this means that a site cannot be overcritical in 
two successive time-steps (for e = the above argument does not work, but the 
statement holds obviously). 

All avalanches start with a single site. Let C\{x) — {i}. Then d(i,j) = 1 for 
all j G C 2 (x). This implies that any two elements of C 2 (i) can be connected 
with a path of length 2, so no two sites of 62(2) are nearest neighbors. If there 
exists a path of even length between two points in A, then all paths connecting 
these points are of even length. Therefore we can repeat the argument proving 
by induction that dA(i,j) G 2Z for all i, j G C n (x). □ 

Proposition 7. The linear maps are all invertible whenever e > 1/2 or 
e > and E c > e/(l — e). // we have E c > e/(l — e), then 

det Lij = e SiJ . 

Proof. Take arbitrary x G R> - First we observe that since the sum over each 
column of Q(x) is less than or equal to 1, we have ||<5(a;)w||i < ||w||i for each 
v G R N . This implies that 

WSixM, = \\(J(x) + (l-e)Q(x))v\U 

> WJWvh -(l-eJHQCrMl! 

> (2e-l)|H|i. 

If e > 1/2, then S(x)v 7^ for all v 7^ 0, so we have invertibility. 

For e = 1/2 the claim follows since in the above chain of inequalities at least 
one is strict if v 7^ 0. In fact, if J(x)v = ev, then Vi = for all relaxed sites 
i. We claim that the equality Q(x)v — v is impossible. To see this denote by 
Q the minor-matrix formed by the rows and columns of Q(x), corresponding to 
exited sites, and denote by v be the respective reduced vector. Then Qv = v. 

Let U be the set of overcritical sites k with Vk = maxu; (we suppose it 
is positive, multiplying by —1 in the opposite case). Choose a boundary site 
fc G U, i.e. the number of neighbors I to fc with vi = Vk is less than 2d. Then: 

Vk = Qkivi < Vk*^2 Q kl - Vk - 
1 1 
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This contradiction yields the result. 

Finally consider the last statement about the case E c > e/(l — e). It is 
proved by reducing the matrix S(x). If x^ < E c , then column Ci(S(x)) equals 
(0, . . . , 0, 1, 0, . . . 0) T , where the 1 is in the i th position. We can start the de- 
composition of det S(x) with column i, and hence we see that row i and column 
i can be removed from S(x) without changing the determinant. We remove 
all rows and columns that correspond to relaxed sites. If p{x) is the number of 
overcritical sites of x, we get a p(x) x p(x) matrix S re d(x). If site k is overcritical 
then Jfcfc(x) = e. If E c > e/(l — e), then it follows from LemmaElthat all nearest 
neighbors of k are relaxed. Hence column k of Q{x) has only zero entries. This 
shows that S rec i{x) = diag(e, . . . , e). Then 

detS(x) =detS red (x) = e p{ - x \ 

so det Lij — e SiJ . □ 

Remark 2. In the original model of Zhang one has e = in which case det Ly — 
if ^ 1. But it is not true that non-trivial kernels can occur for e = only, 
contrary to what was stated in 'BCKjj. A simple counter-example is the case 
N = 2, E c = 1/3 and e = 1/3. For x\ > and 2x\ + 3a;2 < 1 we have: 



Xl 




'2 


3" 




Xl + 1" 


X 2 _ 


)"5 


2 


3 




x 2 



and so det L12 = 0. 

Having non-degenerate maps in the model is more convenient from the point 
of view of mathematical tools (though from a physical viewpoint it can make 
no big difference between degenerate and close-to-degenerate systems). Fortu- 
nately, degenerations occur only for a negligible set of parameters. 

Theorem 8. The maps Lij are invertible for almost all (e,E c ). In fact, they 
are invertible for the parameters complimentary to the set 3 C [0,1) x (0, 00), 
which consists of a finite set of vertical intervals for fixed d and N . 

Proof. Fix an avalanche Ay and let L|- be the corresponding linear maps 
(we stress dependence on e). These maps are the compositions of elementary 
matrices S e {x{r{x, i))) . . . S e (x(l)), with the factors from JTJ 

S*(x(t)) = l + (e-l)(§ -Q)(x(t)) 

being polynomial in e and independent of the choice of x = x(l) G My. The 
condition detL^- = is equivalent to det S e (x(t)) = for some t. Denoting by 
Sp_ (T) the negative part of the spectrum of T, we get: e € 1 + Sp_ (Q — ^f) -1 - 
There are only finite number of possibilities for the matrix S e (x) (though a 
countable number for their compositions L\j , the length of which grow as E c — > 
0). Thus we obtain k = k(d,N) different values of e for which detS* e (a;) = 0: 
{e a }a=i- F° r each e a there is the maximal value E£ of E c (finite if e a 7^ 0), 
where the corresponding matrix S^lx) can appear in the avalanche. Thus the 
set of degenerate systems is {(e, E c ) \ e = e a , < E c < E®}. □ 

By proposition [7| 5 does not intersect the set {e > 1/2} U {E c > e/(l — e)}. 
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2. Removability of singularities and coding 

The map F may be considered as a piecewise affine map F: IxM^IxM, 
where / = [0, 1] and F(t, x) = (Nt mod 1, F [Nt] (x)). The map t ^ Nt mod 1 
is not conjugated to since the points m/N k G / do not have unique repre- 
sentations in T^. However the sets {m/N k } x M C / x M are singularities, and 
following the standard approach for piecewise affine maps, should be removed. 

In some physical systems, like the Bclykh family, the singularities propagate, 
intersecting themselves transversally. The Zhang model is not a general position 
system in this respect, because singularities {m/N k } x M C I x M = M map 
into themselves, forming zero angle. 

2.1. Construction of attractors 

Define the (spatial) singularity set S(F) = UijdMij. Then U = M \ S(F) 
consists of a collection of open connected sets Z = {Z}. Let Uq := U and 

U n :=\jFi(U n -i)r\U. 

ieA 

We say that x € S(F) is a non-essential singularity of order to if there exists 
e > and to > such that 

card{Z ez\u n n B e (x) n Z ^ 0} < 1 

for all n > m. Denote the set of non-essential singularities of order m by 
NES(F; m) and let NES(F) := U m > NES(F; m) be the set of all non-essential 
singularities. Define ES(F) = S(F) \ NES(F) to be the collection of essential 
singularities. Observe that there is a natural extension of F to Vo = U U 
NES(F). In the following we let F denote the extended map. As above we 
define 

V n :=\jFi(V n -i)nVo. 

Let X = (~l n >oK and V = x X. Clearly F(V) = V. The set y = X is called 
the physical (or spatial) attractor of F, and A = x y = V is the extended 
attractor of F. 

Proposition 9. F\x> is continuous. 

Proof. The set T> intersects non-essential singularities only. Hence we must 
show that if x is a non-essential singularity in T>, then the extension of each 
Fi to NES(F) is continuous at the point x. Choose to e N and e > such 
that B e (x) n U n intersects only one partition clement Z G Z for n > to and let 
y G B e / 2 (x) n V. Then B e / 2 {y) C B s (x) and so B e / 2 {y)f\U n intersects the same 
partition element Z. So x and y are mapped by the same affine map Fi\z for 
each ieA. The claim follows. □ 
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In general, the map F does not have a continuous extension to A, but only 
to A\ (S^ x ES(F)). Actually, if x £ ES{F) ny lies on the boundary of several 
continuity partitions for Fj, then there are several extensions of F to (i, x). Thus 
we can continuously extend F to A only when the essential singularities do not 
intersect the attractor (are removable). 

2.2. Symbolic Coding 

If the singularities can affect the dynamics only for a finite number of iterations, 
then the dynamics can be well approximated by a topological Markov chain. 

Definition 2. We say that singularities are removable if there exists m £ N 
such that S{F) = NES(F, m). 

The physically most relevant observables cj> : M — > R arc those that are 
determined by avalanches. We say that (j) is an avalanche observable if it is 
constant on continuity domains [i] x My. 

Theorem 10. If singularities are removable, then the map F is well-defined 
and continuous on A and there is a topological Markov-chain {Y,\,a\) and a 
continuous semi-conjugancy g : A — > such that for all x,y £ A and for all 
avalanche observables (j) we have: 

g(x) = g(y) => 4>(F n {x)) = ^{F n (y)) Vn > . 

The Markov-chain is determined by a matrix A which has a maximal eigenvalue 
equal to N . 

Remark 3. It is clear that all properties related to distribution of avalanche 
size, duration, area, etc. are invariant under a semi-conjugancy such as this. 
Observe that for each avalanche observable 4> on A, there is a unique observable 
<p' : Y>\ — > R such that <p — 4> 9- Suppose we have a measure /i on A, and let 
v = <7*/z. If 4> is an avalanche observable on A, then the statistical properties of 
4> with respect to \i are equivalent to the statistical properties of <fr' with respect to 
v . In this {<j\, S^) is a good approximation to F\j\ L . The coding gives estimates 
on entropy and growth of periodic points, but these estimates are asymptotically 
no better than what we get from the trivial semi-conjugancy M — > S^. 

Proof. Singularities are removable so there exists an integer m £ N such 
that tt s o F m (M) only intersects trivial singularities. Let Xi, . . . ,X S be the 
closure of the connected components of ir s o F m (M). F is well defined and 
continuous on these components. Let Y\,...,Y S be the intersections of the 
components X\, . . . ,X S with y. We construct the partition TZ = {{i]xYk} and 
enumerate it so that TZ = {Ri, ■ ■ ■ , Rr}, where r = Ns. 

Let A = || || be the r x r matrix defined by the rule: ay = 1 if F(Ri)(lRj ^ 
0, and ay = otherwise. A sequence R u)0 R [AJl ... is legal if a Uf _ lUf = 1 for all 
t £ N. Define g : A — > by g(x) = {ujquji . . .ui t . . .), where F t {x) £ R Ut . 
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To prove that g is surjective it suffices to show that for each legal sequence 
Ruj q Ruj 1 • • ■ , there is a point x £ tt s o F m (M) such that F l (x) £ R Ut for all 
i £ N. Note that each oj can be written as a pair (i, k), where t £ {1, . . . , N} 
and fc € {1, . . . , s}. Hence we can write 

oo oo oo 

fl F- n (R Un ) = pi F-([i„] x y fe j = {t} x p F t ; x o . . . o F^(nj . 

n— n— n— 

The continuous image of a connected set is connected, so for each i = 1, . . . , N 
and each k = 1, . . . , s there is a unique I £ {1, . . . , s} such that Fi(Xk) C X;. 
This implies that we have a nested sequence 

y„ c f,; 1 ^) c f,; 1 o F t -\Y kl ) c . . . 

and hence the intersection is non-empty. 

It is clear that gn is continuous (see [J^ for details). Since the partition 
1Z is a refinement of the continuity partition the conjugancy will be injective 
up to the classes of points that follow the same continuity domains. Hence 
if 4>(F n x) 4>(F n y) for some avalanche observable <f> and some n > 0, then 



Remark 4. Suppose we modify the Zhang model by using a full shift (Sat, ct/v) 
as the excitation factor. It is then possible that the modified map F is injective 
on Tin x y . Since we have strict attraction in the spatial factor after a fixed 
number of iterations it is clear that we can then obtain an injective coding, and 
hence a topological conjugancy. However, if we make this modification it is not 
clear that Ejy x y equals the set 



oo 

0= P F n (E N xM). 

n — — oo 

In fact if the maps Fi\y are all injective, then F\q is invertible, but F|j] jv x;y 
is typically non-invertible. The reason for this is that, due to contraction, a 
point x £ y does not have preimages for all the maps F and F ar & invertible 
only on Fi(y) C y. So to obtain invertibility we must turn to the attractor fi. 
From a physical point of view the spatial attractor is of the great interest, so it 
is desirable to have an attractor which is a Cartesian product of the Bernoulli 
shift and the spatial attractor y. 

We can always construct a coding of F\x> (even in non-removable case) by 
choosing a partition TZ = {R\, . . . R r }, and taking gn : T> — * {1, . . . , r} N to be 
the map sending a point x £ T> to the unique sequence u> £ {1, . . . r} N such that 
F l (x) £ Ru t for all t > 0. But there is no reason, however, to expect gn(P) to 
be a topological Markov chain, cf . |BCK| . 
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3. Metric properties 

The natural volume on M is given by the product measure of the uniform 
Bernoulli measure on E^ and the Lebesgue measure on M. By iterating this 
measure (and averaging) we can construct SRB-measures. However it can hap- 
pen that the measures constructed are supported on essential singularities, 
where it is not possible to define the dynamics in such a way that the mea- 
sure is invariant. Hence we must give some conditions to ensure the existence of 
SRB-measures. If there is an SRB-measure it is characterized by the fact that 
its projection to E^- coincides with the uniform Bernoulli measure. From this 
it follows that any SRB-measure is a measure of maximal entropy. In situations 
where the system allows symbolic coding the SRB-measure corresponds to the 
Perry measure on the topological Markov chain E ^ . 

3.1. Existence and characterization of SRB-measures 

Let m = m u x m s , where m u = /i Bor is the uniform Bernoulli measure on E^-, 
and m s = /i Lob is the Lebesgue measure on M. We say that an invariant Borel 
probability measure /ionM has the SRB-property if there exists a measurable 
invariant set G C M such that 

1. m{G) > 

2. m u {v u (G)) = 1 

3. All points x £ G are future generic with respect to [i, i.e. 

n-1 

for all x e G and all continuous functions <j> : M — ► K. 

For Axiom A attractors one can ensure the existence of measures for which 
the set of generic points has full Lebesgue measure and this is equivalent to 
saying that the canonical family of conditional measures on the unstable man- 
ifolds are absolutely continuous with respect to the Lebesgue measure. For 
non-invertible maps one can in general only expect the set of generic points to 
have positive measure and hence it is unreasonable to require that m{G) = 1. 
Condition 2 is (for physical reasons) important in the Zhang model. It means 
that the statistical properties do not depend on the choice of a generic sequence 
t of excitations. (This is always implicitly assumed in the numerical investiga- 
tions of the Zhang model that can be found in the physical literature.) Moreover, 
condition 2 will be satisfied for the SRB-measures that can be constructed by 
iterating the measure m. 

By a standard approach we can give conditions for existence of SRB-measures 
that hold if singularities are removable, but it is not known if these conditions 
hold in all non-removable situations. 
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Proposition 11. Let (e,E c ) does not belong to the negligible set S of Theorem 
[3 If there exists n > 0, C > and q > such that 

V6 > 0,Vt > : m(F-'(S+ x U s (ES(F;n))) < CS q , 

then there exists a set T> C M ( constructed in ^2. 1)) . which may intersect singu- 
larities, and a natural extension of F to T> such that F(T>) = T>. Moreover the 
set T> carries an F -invariant Borel probability measure with the SRB-property. 

Remark 5. Provosition M 1\ is a simple modification of the result of Schmeling 
and Troubetzkoy JfSTf . In their paper the conditions for existence are in general 
too restrictive for the Zhang model. In fact, in Example A of ^ we show a 
situation where the SRB-measure constructed in \STj does not exist, but we 
clearly have existence of a physically relevant measure. The reason for this 
paradox is that one in general remove all singular points on the construction of 
the attractor, even if there is a natural extension of F to the points of singularity. 
( Provosition M 1\ obviously applies to this example since ES(F;5) = 0J 

Proof. In |ST| it is shown that a piecewise smooth map / with singularity set 5* 
has a measure, not supported on singularities, such that the set of generic points 
has positive Lebesgue measure. They require that the following conditions are 
satisfied: 

1. The restrictions of / to each of its continuity domains are diffeomorphisms 
onto their image. 

2. The second differentials D 2 f x does not grow too fast close to singularities. 
(See |ST) for a more precise formulation.) 

3. / is hyperbolic. In this context this means that there are constants C > 
and A € (0, 1) such that for all x $ S there is a splitting of the tangent 
space at x into subspaces E + {x) and E~(x). There are cones C + (x) and 
C~ (x) around E + {x) and E~(x) that are invariant under Df x and Df^ 1 
respectively. The angles between the C + (x) and C~(x) are bounded away 
from zero, and for all points x that do not intersect singularities in the 
first n iterations it holds: 

||-D X /»|| > C- 1 A-*»]]t7|| for v e C+(x) , 

and 

||-D x / n (w)|| < CA n |M| for v e C~(x) . 

4. There exists C > and q > such that m(/-*(f7 e (S))) < Ce q for all 
e > and all t € N. 

We apply this result to the map F\^+ X ^ u ^s(F-n))' singularity set for 

this map is contained in S^- x ES(F), so by assumption condition 4 is satisfied. 
Condition 1 follows from Proposition0 condition two is obviously satisfied since 
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F is piecewise affine and condition 3 follows from Theorem[5]with E + = M 1 ©0, 
E~ = © 1* and C ± being the regular cones around them (actually Theorem 
[S] ensures hyperbolicity for some iterate F T , which implies the claim). 

In |ST| the measures are constructed by iterating m, averaging and taking 
a weak limit. It is clear that, in the Zhang model, any measure obtained in this 
way will satisfy condition 2 in our definition of an SRB-measure. □ 

If an SRB-measure exists it can be characterized by a number of different 
properties. From a physical perspective it is reasonable to require that a relevant 
invariant measure should preserve the uniform Bernoulli structure on E^- . This 
corresponds to the Lebesgue measure on [0, 1] in the alternative formulation of 
the map F, and hence to absolutely continuous measure conditional measures 
on the unstable space [0,1]. 

Proposition 12. If fi is an SRB-measure on T>, then fi u :— (7T u )*/i is the 
uniform Bernoulli measure on E^. 

Proof. There is a set A — tt u (G) of full m u -measure, such that all t 6 A are 
generic with respect to fi u . Take a continuous function <p : E^- — > K. Then 

/ 6du v = lim -V#+)'t)= / c/)dm u , 

J ri^oo n ^ J 

where the left equality holds for t d A and the right one for t S B with B C E^ 
a subset of full m^-measure (from Birkhoff ergodic theorem). Since AC] B ^0, 
we get: J cf>dfi u = J (f>dm u for all continuous functions <p. □ 



3.2. Measures of maximal entropy 

Suppose that there exists an invariant Borel probability measure \x on M . Let 
[i u := (~k u )*H and let {ft} to be the canonical family of conditional measures 
on the fibers ^"""^({t}). By the Abramov-Rokhlin formula 

h^(F) = h^(a+) + h,(F\a+), 

where 

K{F\a + N ;Q)= lim - / H„ t ( \/ o • • ■ o F^-^Q)) d^(t) , 

n— >oo n / V v / 

for a partition Q and 

h IM (F\a+)= sup h^F\a+;Q), 
Q 

The supremum is taken over all finite measurable partitions Q of M. This 
formula was originally proved for product measures by Abramov and Rokhlin 
.AH! and extended to arbitrary skew products by Bogenschultz and Crauel BC . 

Below we use the notation _Ft for the dynamics over a pre-fixed sequence 
t = (toti . . . ) G E^y. By n-the iteration we mean the map Fp = F tn _ 1 o • • • o F to . 
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Theorem 13. Iffi is an invariant Borel probability measure on M , then h fl (F) = 
V(cr+). 

Proof. We prove the proposition by estimating h fJ ,(F\a'^ ! ) from above. Let Q 
be a partition of M and 

n— l n—1 

Qh...tn-r ■= V (^-i ° • • • ° ^o) _1 (C) - V ^"'(Q) • 
fc=0 k=0 

Fix e > and choose the partition Q such that h fJi (F\a~^; Q) + e > /^(i^cr^). 
The maps Fj are non-expanding so it follows from the Ruelle-Margulis inequality 
[EH] that 

-H vt [ \f F- k (Q)) -^ hut (F t ) = 0. 
k=a 

Moreover the convergence is //"-uniform and so the same holds for the integrals. 
Another way to see it is via the multiplicity notion of i |4.2l (then Q should be 
subordinate to each continuity partition {M^ \j = l,..., qi}): 

hu(F\at;Q) < lim - logmaxmult(Q to ... t „_ 1 n supp(i/ t )). 

n— >oc n |t|<n 

Therefore h f _ l (F\cr+) < e. Let e -> 0. □ 



Remark 6. Theorem M'A is a partial case of Theorem 4 from IKR2I. 

We say that an invariant measure fi is maximal if h^F) — sup„ h v (F), 
where the supremum is taken over all invariant Borel probability measures on 
M. It follows from Theorem US that h M (F) < logN. So if h top {F) > log TV, 
the variational principle fails (this can happen for piece- wise afhne systems, see 
|KR1I l?vR2p . But we show in ^4.21 that the abnormal growth of ht op (F) does 
not occur in the Zhang model, at least for generic values of parameters E c , e. 

Corollary 1. Any SRB-measure on T> has entropy h^(F) = log N and is hence 
a maximal measure. 

Corollary 2. Suppose singularities are removable and that fj, is an SRB-measure 
on A. Let g : A ~> Ej^ be the semi-conjugacy constructed in the proof of Theorem 
\1(A If (<xt, Ej[) is topologically transitive, then g*/i is the Perry measure on Ej[. 

Proof. A transitive topological Markov chain has a unique measure of maximal 
entropy. This measure is called the Perry measure jKHj . □ 
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3.3. Hyperbolic structure 

There are several ways to define Lyapunov exponents for the Zhang model. The 
Zhang model can be represented as a piecewise affine map, where Bernoulli shift 
is represented as the expanding map t i— > Nt mod 1 of the interval (see H4.2fl . 
Hence it is clear that there is one positive Lyapunov exponent Xo = ^ogN. We 
define the other exponents by introducing the co-cycle T : M — > GL(N,M), 
defined by T(x) = L ij: where x £ [i] x M tj . For x £ M and u G M. N \ {0} we 
define 

X{x,v) = Bin l^mga^WM, 

ra— >oo n ||u|| 

It is a general fact that the function x(i, ■) takes at most N different values 
Xi(x)>---> Xn{x)- 

Proposition 14. For all x £ M £/ie Lyapunov spectrum is: 

xt = log AT > > xr( £ ) > • • > Xw(^) 

and for (e, E c ) outside the negligible set S /rom Theorem^; Xjv(^) > — oo. 

Proof. From Theorem [3] we know that there exists T £ N and c € (0, 1) 
such that 

||r(F T - 1 (x))...r(^(x))r(£)||< c 

for all x £ M. It immediately follows that xipi v ) — T^ 1 logc < 0. 

For e > 1/2 and arbitrary E c or for e > and i? c > (1 + e)/(l — e) all 
linear maps are invertible, and so for all x £ M and all v £ R N \ {0} we have 
x(x, v) > logfc, where fc = min^ minSpfXjj) > 0. □ 

If there exists a unique SRB-measure, then it follows from the Osceledec 
theorem that there are numbers Xi i ■ ■ ■ > Xn such that Xi(%) = Xi f° r Lebesgue 
almost every x £ M. The numbers Xo i Xii ■ • ■ i Xn are the Lyapunov exponents 
of the Zhang model. From Proposition [21 it follows that the Zhang model is 
hyperbolic in the sense that the Lyapunov spectrum consists of: 

X([ = log TV > > Xi > ■ ■ ■ > X N > -oo . 

We see from Corollary ^ that the Pesin formula h fJi (F) = x + holds for any SRB- 
measure. If there is no SRB-measure then the Lyapunov spectrum should be 
defined as functions on M: 

Xi{x) = / X»( t I 2 ; )^Bo r , 

where /i B cr is the uniform Bernoulli measure on E^. 
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3.4. Entropy of physical vs. mathematical models 

We can reformulate the Zhang system as the map / : B — ► B, where B = 
l ->i>of l (M) is a compact /-invariant subset of S^- x M. N . We wish to compare 
this to the induced transformation F : M — > M (cf. Remark 0. 

If /i is an F- invariant Borel probability measure on M, then there is an associ- 
ated /-invariant Borel probability measure fi on B (and vice versa). Abramov's 
theorem (\Br\) relates the entropies of both systems: 

hfi(f) = h^F) ■ A(M). (3) 

One does not need to assume ergodicity and can allow degenerations |DGS| . as 
happens for the case of Zhang model. In ergodic situation by the recurrence 
theorem of Kac |Brj for a /i-generic point x £ M: 

1 1 ™ _1 

„_= ii m -X" T (F k x), (4) 

A(M) n ^ 

where t(x) — r(i, x) is the avalanche time initiated by addition of e, to x G M 
(see ^1.2|) . If n = /isRB is a unique SRB-measure, the above point x can be 
chosen Lebesgue generic. The resulting limit is the average avalanche time f 
(we discuss it in more details in tl7.1l7.2|) and we obtain: 

%(/) = h^(F)/f resp. hp, SRB (f) = h ^ / ? . 

In general non-ergodic situation to get equality we should integrate the 
terms in right-hand side and then we again obtain the average avalanche size 
(r), but now it is the space-average. Substituting this into © we get: 

h fi (f) = h fi (F)/(r). (5) 

For SRB-measures this formula is indicated by the Ledrappier- Young theo- 
rem |LY| . because we have only one positive Lyapunov exponent. 



4. Topological entropy 

To calculate the topological entropy of F we established in |KR2| a set of in- 
equalities, using the technique developed by J. Buzzi |B1| . |B2| for piecewise 
expanding maps and piecewise isometries, see Appendix^ The contraction in 
the maps F^j provides difficulties, so several results were generalized to fit the 
framework of the Zhang model. It is not however true (as was widely believed, 
see |B2| ) that the contraction does not contribute to topological (contrary to 
metric) entropy, the corresponding counter-example can be found in |KR1| . The 
Zhang model has a feature common to all such examples KR2 , namely angular 
expansion, but still for most values of the parameters this abnormal increase of 
the entropy does not occur. 
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4.1. Growth of the number of continuity domains 

Let V = {[i] x Mij} be the partition of continuity for F, and enumerate the 
elements so that V = {Pi, . . . P r }. Let 

n-1 

[Poo ■■ -Pan-,] ■■= f| F ~ m ( P "m), 

m=0 

and 

V n = {[P ao ...P an _ 1 ]^®\a m = l,...,r}. 
We define the singularity entropy of F by 

H sing (F) = lim -logcard(P"). 



Remark 7. Define a map g : A — > by letting g{x) be the unique sequence 
a ai . . . such that F n (x) e P an . Then 

H sing (F) = h top (a+;g(A)). 

For piecewise affine expanding maps and piecewise isometries it is clear that 
this also equals the topological entropy, but due to the contraction, this is not 
obvious in the Zhang model. In addition, if singularities are not removable, the 
map g has discontinuities. 

Definition 3. Call a point x e S{F) an unstable singularity if for all i and 
k ^ I we have: \im y ^ x F ik (y) ^ \im y ^ x F a (y). 

Theorem 15. // all singularities S(F)C\y are unstable, then h top {F) — H s i Dg (F). 

We need the following technical lemma: 

Lemma 16. // the singularities in y are unstable, then there exists a constant 
7 > such that for all 5 > 0, x e y n M ik and y e y n Mu, k ^ I, we have: 

d(x,y)<S^d(F i (x),F i (y))>6. 

Proof. Suppose that for all 5 > there exists x € M ik n y and y € Mu n y 
such that d(x,y) < 5 and d(Fik(x), Fu(y)) < 8. There exist sequences {x m } C 
Mik and {y m } C Mu such that d(x m ,y m ) — > and d(Fik(x m ), Fu(y m )) — > 0. 
The sequence {x m } has a convergent subsequence x mn — * z. The point z lies 
in S(F) and y mn — > z. By the continuity of the maps Fik and Fu we have 
Fik(x m „) — * Fik(z) and Fu(y m „) — * Fu{z). Since the metric d is continuous on 
M x M we have 

d(Fifc(a:), FaCO) = lim ^aK),fii(K)) = 0- 
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Hence F^(z) = Fu(z). Contradiction. □ 
Proof of Theorem 1151 Let 

[P ao . . . P 0n _J = [i ...in-i] x K, 

where K C M is a convex polygon. Fix <5 > and set k = [log 1/5]. Let 
Z\, . . . , Z m r$\ be a (5-spanning set for K. Chose t £ and define ^V fe sequences 

s r ...r fc _! = (io • • • ira-lXo • ■ ■ ^fc-l^n+fcin+fc-l • • ■ ) S . 

Since the maps are contracting it is clear that the set {x = (s ro ... r)k _ 1 , 2;)} 
is a (n, 5)-spanning set for [P ao . . . P an _ 1 ], and since the minimum number of 
balls needed to cover a convex polygon K C M is bounded by m(8) < Cn5~ n 
we see that the number of (n, <5)-balls to cover [P ao . . .P a „_ 1 ] is bounded by 
m(S)N k ■ #{[-P Qo . . .-P a „_i]}- Therefore we get an estimate for the number of 
(n, (5)-balls to cover M and so 

log(C^^card(P»)) 
htop(F) < km i '- = H sing (F) . 

n— *oo 7| 

To see the opposite inequality let A,B£ V n and take x\ = (t, x) £ A and 
x 2 = (s, y) G B. Suppose that A ^ B and that to = s 0i • ■ • i tn-i — s n-i- Then 
there is m < n such that n s o F rn (x) £ M tm ^ and 7r s o F m (y) £ M tm i, k ^ I. By 
Lemma [TBI there is 7 > such that for all £ < 7: 

max{d(7r s o F m {x),* s o F m (y)), d(n s o F m+1 (i),^ ° P m+1 (y))} > € • 

Therefore for <5 sufficiently small, no (n + 1, <5)-ball can contain points of both ^4 
and £? and the minimal (n + 1, <5)-spanning set has at least card(P™) elements. 
Then h top {F) > H sing (F). □ 

Theorem 17. For the Zhang models: h top (f) = H s i ng (f), h top (F) = H s i ng (F). 

Proof. By Remarkdall quantities are topological entropies. The corresponding 
systems in the second equality are the Poincare return maps for the transfor- 
mations of the first equality. The map F : M — * M is already the return map 
for / : B — ► B by the very construction. 

To achieve the same claim for the symbolic system we extend the partition 
V of M to a partition V of B, which on £+ x (R^ \ M) n B equals the product 

of the standard partition = L>i[i] and the partition of the spacial part by the 
hyperplanes {xi = E c }. Denote by s < r + N(2 N — 1) the number of elements 
of the new partition V '. 

Let B C B be the /-invariant closure of A in x R> . Define a map g : 

B — ► by letting be the unique sequence 6061 . . . such that f n (x) £ Vb n - 
We wish to prove that 

h top (f) = H sing (f) = h top (crf\g(B)) ■ (6) 



22 



Dynamics and entropy of SOC 



For this it is sufficient to check that the singularities S(f) fl tt s {B) are unstable 
for / (the second equality follows from the definition). 

pN 



Consider a singular point x £ K> . Let y,z tend to x by two differ- 



ent domains of the projected partition V in K> . Then we can subdivide 
{1, . . . , N} = A U B U C, where y t = E c + 0, Zi = E c - (with the obvi- 
ous notations instead of limits) for % S A, y.- L — E c — 0, Zi — E c + for i e 5 
and Zi belong to the same side of -E c for i G C. Then for 5 = — S'(.z) we 
have: 

e-l, if i = j € A, 
1 - e, if i = j e S, 



Sij — { 



^f, ifd A (i,i) = l,j€A 
i{d A (i,j) = lJeB, 
otherwise. 



We should check that S ■ v ^ for any vector v £ M> with components Vi = E c 
for i E A U B (this implies that a; is unstable). Note that other components Mj, 
i £ C do not contribute to the product. 

Let i be a site from A with not all neighbors from A (if the set A is empty 
consider B). Denote the number of yl-neighbors of i by fc^ < 2d and the number 
of B-neighbors by k B > 0. Then (S ■ v), = E c (2d - k A + k B )^± ^ 0. 

Thus Theorem 1151 implies (jHJl. Moreover, the same reasons yield a more 
general statement. Namely, since g is a semi-conjugancy, we get: 

ho P (f\}Q=h top (at\9(fC)) (7) 

for any subset KcB (recall that g may have discontinuities, but our arguments 
are not injured by this fact). This subset needs not to be invariant, and in this 
case we should use Bowen's definition of entropy |Boj . Thus Katok's entropy 
formula ^ implies that hp,(f) = hg t p,(a^) for all /-invariant measures jl, which 
are not supported on singularities. 

If the system (M,F) possesses a measure /i of maximal entropy, then by 
the obtained result the second claim of the theorem follows from © and Kac's 
theorem |PY| . In general, we can apply the above arguments to the partition 
of M by the subsets of equal return times and using the fact, that both returns 
of / and g have the same combinatorics, we get: h top (F) = h top (u^\g(A)) (see, 
for instance, the loop equation approach |Ptp. □ 



4.2. Evaluation of topological entropy 

It was predicted on the base of variational principle in jBCK] that topological 
entropy of the Zhang model is htop(F) = log N. However this principle does not 
apply because the map is not well-defined (continuously) on the whole space (or 
thanks to non-compactness if we remove the singularities). In fact, there can 
be no invariant measures on the non-singular part at all. 

While we support the claim that h top (F) = logA^, it will not be proved in 
full generality. We start with the asymptotic statement. 
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Figure 1: Shows the avalanche type domains on (e, i? c )-bifurcation diagram for N — 2. 
The top domain is E c > j^- . The line from infinity to the origin is E c — rr^— . We see 
infinitely many domains with different avalanches, accumulating to two of the axes. 



Consider the bifurcation diagram on the (E, e) strip (0, +oo) x [0, 1), where 
a point is critical if in its neighborhood dynamics of the Zhang model can expe- 
rience avalanches of different types. Thus the strip is partitioned into different 
avalanche type domains. The partition depends on L, d. For N — 2 the diagram 
is shown on Figure ^ 

Note that for all L, d there is the top avalanche type domain representing 
the shortest avalanche time. For N = 2 it is given by the relation E > j^-. 
Also note that some domains have e-projection strictly smaller than the interval 
[0, 1). The next statement concerns only the top domain and the avalanche type 
domains that are adjacent to the line e = 1. 

Theorem 18. For the Zhang model: < h top (F) — logiV < 9(E,e), where 
liniE^oo 9(E C , e) = (e fixed) and lim e _»i 9(E C , e) — 0. In the latter case E c 
changes accordingly with e so that (E c ,e) belongs to the same avalanche type 
domain (then E c — > oo, though this case differs from the former). 

We can remove the TV-rational points from and then represent it as 
the subset of / = [0,1] with points n/N k being deleted (this process does not 
change the topological entropy). In fact, we need to remove only points n/N 
for the others will be deleted by inverse iterations of the map oi I x X with the 
formula F(t,x) = (Nt, fwt] ( x ))- Thus we represent our system as a piece-wise 
affinc partially hyperbolic map of the subset of R 1+Ar . This is important for an 
application of the results from Appendix lAl 

The proof uses the notions of multiplicity and multiplicity entropy, due to 
G. Keller and J. Buzzi. Given a finite partition Z = {Z^} of M we define 

z n = {[z ...Z n - 1 ]^<D}, 
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where 

[z ... z n _i] = z Q n F~\Zx) n • • • n F^iz^x) . 

Define multiplicity of a partition by 

mult (Z) = maxcard{Z e Z \ Z 3 x} , 

If Z is the continuity partition of the map F we often denote the multiplicity of 
the Z by mult(F). Then it is clear that mult(F") = mult(Z"). The multiplicity 
entropy of F is (the limit exists by subadditivity, cf. KH ) 

H muit (F) = sup lim - log mult (Z n ) , 

and we see that if Z is the continuity partition, then 

H n -MF) = lim -logmult^™). 

n — >oo 77, 

It is clear that H mu \ t (F) = if the singularities are removable, so h top {F) = 
\ogN by the result in Appendix^] For big E c 3> 1 the singularities are gener- 
ally non-removable, but still we have the same effect asymptotically: 

Proof of Theorem 1181 We take 6(E C , e) — H mu \ t (F). Since the singularities 
t = n/N £ I of the map er^ : 1 1— » Nt do not intersect in inverse iterations, the 
multiplicity growth is only due to the spacial maps Fi : M — > M . Thus using 
the notation F t from Section 13.21 for the dynamics over a prefixed sequence of 
excitations we obtain H mu \ t (F) = sup t6S + i? mu it(^t)- 

To show the first claim let us notice that when e = const, but E c — > oo, then 
the avalanches map the critical part of the boundary {3i : Xi = i? c } C dM far 
from dM, namely to the distance ~ 7(5, L, d)E c , see Figures [TBI for N — 2 and 
1 171 for N = 3. To reach again the boundary and experience avalanche we need 
many shifts by the basic vectors ej. Thus the singularity can meet only after 
big number of iterations. Since the initial picture of singularities has bounded 
multiplicity (for (e,E c ) from the top avalanche type domain), the multiplicity 
decreases at least as k/E c so that it vanishes in the limit. 

To prove the second claim we use the inequality H s i ng (F) < J2Pi(F) h° m 
Appendix IA.3I If the avalanche type domain is fixed, the number of compo- 
sitions of matrices S(x) in one avalanche (see HI. 3(1 is bounded. Every such a 
matrix tends to identity when e — > 1 . Thus all the linear parts of avalanche 
maps Fij tend to identity and so angular expansions Pi{F) tend to zero. □ 

Notice that if E c is fixed, but e — * 1, then the number of avalanches has 
unlimited grow and the previous argument do not work. However, due to esti- 
mates of i jl.2l on the maximal avalanche length r m we conclude that 6(E C , e) — > 
as either E c — * 00 or e — ► 1, both quantities being related by the constraint 
E c > N(l - e)~? diam(A)-cr £ Qr somc a > q ( we don't require, but allow TV — * 00 
as well). This statement is stronger than in the theorem. 

Now we are going to prove vanishing of 9(E C , e) a.e. in the finite part. 
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Theorem 19. For generic (E c ,e) we have: h top (F) — log AT. 

We will prove the statement not only for the Zhang model, but also for nearly 
Zhang models. By this we mean the following. The map F is a bundle over the 
Bernoulli shifts <t^ with factors Fi being piece-wise affine partially contracting 
maps. That is M — LLfMy and for = Fi\M i:j we have Fij(x) = Lij(x) + bij, 
bij = Lijei. We are going to make arbitrary small generic perturbations of the 
matrices Lij (still with spectrum within unit ball) and vectors bij (one should 
care that My are mapped into M) and prove the statement for this modified 
system. 

Due to round-off errors there is no much difference between the original and 
the perturbed systems in computer simulations. And from the point of view of 
the experiment such perturbations (instrument instability) are indispensable - 
look at |Ru| for the discussion of physical relevance of variation of parameters 
as the noise. 

Proof of Theorem 1191 We claim that for generic (E c ,e) the singularities do 
not multiply. Actually, some intersections of singularities are deformable as we 
vary the parameters, just by the transversality reasons, but the other disappear 
with small perturbations. 

Namely, for a multi-index a = (a±, . . . ,<Xk), ol s = (i s ,j s ) coding an orbit, 
denote F a — F ak o- • -oF ai the corresponding map along the orbit. The singular- 
ities of this map are: Sing(F CT ) = U^ =1 F~^ Sing(F ar ), where o [ r ] = (a\, . . . , a r ). 

If zq € Sing(i 7 ' (T ), then its orbit (with multi-possibilities due to singularities: 
mapping a singular point we extend the components of the map F a in various 
ways) meets several singularity planes, i.e. for some cuts a s — <7r r 1 of a, s = 
1, . . . , m, we have the following system: 

F ai {z ) = zi, . . . , F am (z ) = zi, l qi (zi) =0, lg m (z m ) = 0, (8) 

where l qs (z) are equations for the singularity hyperplanes of the corresponding 
map Fi (there are also inequalities, which we don't mention). When m < N 
there are occasions, when (jSJ has a solution continuously depending on (E c , e). 
However for m > N this is no longer the case. In fact, considering nearly-Zhang 
models we see that for generic data the above system (JHJ is characterized by a 
collection of non-trivial polynomial in e equations (for each E c ). 

More precisely, the set of (iy,6«) giving trivial polynomials has positive 
codimension and hence zero Lebesgue measure. Uniting these sets over all 
choices of multi- indices (ui, . . . , o~ m ) we see that the complement has full mea- 
sure and dimension and so a generic perturbation yields the data (Lij, bij) from 
it. Since non-trivial polynomials have only finite number of zeros, then for 
a generic nearly Zhang model and every E c there is a countable subset of 
{0 < e < 1}, so that the corresponding systems (JHJ) have no solutions. This 
means that multiplicity of F^ does not grow with n and so H mu \ t (Ft) = 0. 

Now let us look to the Zhang model. We restrict for simplicity of exposition 
to the case of two sites A^ = 2. In this case the linear parts of the affine maps 
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are compositions of L\ = 1 + (e — l)Ai and L 2 = 1 + (e — 1)^2 with 



Ax = 



1 

-1/2 



A 2 = J~ l A x J = 



-1/2 
1 



where J 



1 



(we exclude the obvious matrix 1, see Example A for details). Notice that 
det L\ — det L 2 = e. 

Suppose that JHJl has a continuous solution in some domain of (E c , e). Then 
it is algebraic in e and linear in E c (the latter is because the singularity lines 
within one avalanche type domain shift with velocity 1 in the direction of either 
'1\ /V 



ei = 



or e 2 = 



. Differentiating this by E c we obtain: 
Lpi ( z o) = z \i ■ ■ ■ > L Pm (z ) = z m , 



(9) 



where L Pi — L Pi t . o • • • o L Pi 1 for pi — (pi,!, . . . , Pi,ti) is the linear part of F ai 
(pi is different from <Xj because is a composition of several maps L s ). The 
points Zfc, 1 < k < to, are constrained to the singularity lines and we can suppose 
these are the lines {x\ = E c } or {X2 = E c } (all other singularities are mapped 
to them within one avalanche). Thus z' k = Vk + V'feWfe? where Vk = e± or ei and 
Wi = Jvi is the other basic vector, ipk being an unknown scalar. 

Let £fc = L~^(vk), Ck = Lp^iwk); these vectors depend meromorphically on 
e. System © is solvable iff the affine lines + tpkCk, 1 < k < m, in M. 2 have a 
common point. We can suppose to = 3. 

Denote by 0(£, 77) = ( J£, 77) the standard symplectic form on M 2 . The above 
3 lines intersect jointly iff 



0(6,CiMC2,C3) 

t3 



n(6,Cj)n(C3,Ci) 



Dividing by rife=i^(?fe'Cfe) an d using the fact that fi(£fc,Cfc) 
£ -|pfcl we g e t the equivalent equation: 



n(f3,C3)fi(Ci,Ca) = o. 

det(i^) 



^(771,772) + 0(772,773) + ^(773,771) = 0. 
Here 77^ = e^^ Pk ^L^Wk = L Pk Wk, where Z r = L Tl o • • • o L T 



for t = 



(10) 
.,r t ) 



and Lfe = e L fc = 1 + (e — l)Afe is the adjunct matrix for Lk, which gives 




1/2 



1/2 




Now equation (|10fl holds iff there exist 3 not simultaneously zero numbers 
lh>fh>lh such that 0i + (5 2 + P3 = and ft 771 + /3 2 772 + /3 3 77 3 = 0. Since 77^ is a 
polynomial matrix of degree \pk \ and the products of A t are always proportional 
to e\ or e 2 (depending on the left-most factor) this last equation is never satisfied 
if the multi- indices pk are different. □ 

Remark 8. To support the usage of nearly-Zhang models note that the whole 
paradigm of SOC should allow generic perturbations of the data, for if there is 
a fine tuning of parameters, the model is unappropriate for physical explanation 
( of course, we should pass to the thermodynamic limit, but in practice this only 
means some large finite parameters) . 
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Our computer experiments did not expose any exponential growth of multi- 
plicity in the Zhang model (though we see growth in complications of singular- 
ities), so we suggest that H mu \ t (F) — and hence h top (F) = log TV always. In 
addition, by the above discussion we can disregard these exceptional values of 
(E c , e) even if there are any. This finishes discussion of topological entropy. 

5. Geometry of the attractor 

The construction of the spacial attractor y can be interpreted as an iterated 
function system (IFS), where the maps Fi are not affine as usually considered, 
but piecewise affine. Hence one might expect that attractors y are fractal, 
but with various size characteristics, like dimension and measure, depending on 
parameters E c and e. 

5.1. Fractal structure 

Computer experiments show that in certain cases the spacial attractor y has 
fractal structure, see e.g. Figures and |3 We clearly see that y consists of 
self-similar pieces. However the pieces overlap, making evaluation of the fractal 
dimension difficult. So we can provide only estimates of the attractor's size. 

Nevertheless we observe from our experiments that Hausdorff dimension 
dimn(3^) and the Lebesgue measure /LtLeb(3^) of the attractor grow piece- wise 
monotonically with E c and e. Thus the following effects occur in steps: 

• The dimension and the measure of y vanish. 

• dining) is positive, while the measure is zero. 

• Both dining) and Hheb(y) are positive. 

• The attractor y contains an interior point. 

We will demonstrate the dimensional part in the next section, while we 
disregard the observation about the measure. The reason for this is that /iLcb 
is not physically motivated and we should look for an SRB-measure. 

The experiments show that such a measure exists and has support lying 
strictly inside y. (See Figure 0] and Figure 03 In Figure 0] the attractor y is 
shown for the case N — 2, E c — 20 and e = 2/3, and Figure shows the orbit 
of a random initial condition. The latter corresponds to the support of the 
SRB-measure.) This is possible because the contraction rate of ft is smaller for 
the exceptional sequences t € S^, than for a generic one. Thus study of the 
IFS-attractor does not lead to conclusions about ergodicity or uniqueness of the 
SRB-measure. Still it provides an information about spacial distribution of the 
orbits in the Zhang dynamics. 
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Figure 2: Shows the set U w for N = 2, E c = 5 and e = 1/2. 



/ w 
/ \/ 









Figure 3: Shows the set U w for N = 2, £ c = 3 and e = 1/5. 
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Figure 4: Shows the set y in the top right corner of M for N = 2, E c — 20 and 
e = 2/3. The dimension of y is 2 in this example. 



5.2. Dimensional study of the attractor 

The fractal properties of y do not hold for all values of parameters (E c , e). An 
example where y has integer dimension is shown in Figure 01 

It was noted in jBCKj that Hausdorff (fractal) dimension of the attractor is 
about to increase as E c grows. The arguments were the following: For bigger E c 
the contraction rate decreases, so the theory of iterated function system (IFS) 
implies increasing of the Hausdorff dimension 

®y[e,E c ) =dim H (y) 

as a function of E c . While this seems to be true, the statement does not hold 
in precise sense. For instance, for N = 2, E c g [jit, jzr^] the attractor is the 
set of 3 points, while it seems to have non-zero dimension for other parameters 
(computer simulations clearly show this). 

The problem is with the framework of IFS, where usually only conformal 
maps are considered and certain regularity of their mapping graph and overlaps 
is assumed. However we will show validity of the claim in the asymptotic sense: 

Theorem 20. With fixed d, L and generic e we have: lime,.— k» Dy(e, E c ) = N. 
Moreover, T>y(e, E c ) — N for big values E c 3> 1/e. 

On the other hand for all e € [0, 1) it holds: limB c ^o (e, E c ) = 0. 

In the above statement "generic" means both full Lebesgue measure and 
second Baire category. In fact, the equality holds for all e outside a countable 
set. It seems though that the limit statement is valid for all e. Thus we see that 
T)y(e, E c ) is not strictly monotone in E c for its large values as one might expect 
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Figure 5: Shows an orbit of a randomly chosen point (i.e. the support of the SRB- 
measure) for N = 2, E c = 13 and e = 2/3. 



from the arguments cited before the theorem. 

Proof. Consider first the statement about big energies E c ^> 1/e. Let us start 
by demonstrating the idea of the proof on the example N = 2, see Figure IT^l 

The image of the vertical continuity domain M13 of height 1 adjusted to the 
right-top corner is a trapezium with the slope depending on e (see <|11[) of Section 
|6l for the numeration of domains). It is thin - of constant length of horizontal 
section equal e near its bottom side, but long - with diameter approximately 
equal E c l ~^ . Thus if we shift ~ C/e times this domain up and then all of the 
shifted images horizontally to the right, so that its first coordinate satisfies the 
inequality E c — 1 < x\ < E c , then these shifts cover an open domain, including 
a unit square, in the vertical strip K\ — {x\ £ (E c — 1,E C ]} (note that we can 
leave a copy of the domain since this corresponds to the shift - dropping of 
energy (0:1,3:2) | — > {xi,X2 + 1) on M13 C K\ before the avalanche). An easy 
calculation shows that such shifts cover the whole upper part of K\, strictly 
including the continuity domain M13 C K\ adjusted to (E C ,E C ) and covering 
a vertical part of M 12 . A similar scenario happen to the second coordinate. 
Thus in iterating the dynamics we will always have two continuity domains M 13 
and M23 adjusted to (E C ,E C ) and the adjacent parts of M 12 ,M 2 2 lying in the 
attractor. 

For general L and d we observe the same picture: With generic e a continuity 
domain Mu adjusted to the upper-most corner is mapped under avalanche-map 
Fij to the trapezoid-like polyhedron with irrational slopes. Its shifts cover then 
an open domain in each of the strips K t = {x.i £ (E c — 1, E c ]} and so after more 
shifts - the upper part of this strip, whence the statement. 

We illustrate this process on Figure El The 3 domains adjusted to the 
corner (E c , E c , E c ) are mapped into interior of the cube and they have different 
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irrational slopes (we picture them of zero thickness that corresponds to large 
values of E c ^> 1), so that their shifts cover a big open domain near the faces 
adjusted to the above corner. 

Consider now the second statement, £ c <Cl, To estimate the fractal dimen- 
sion from above we use the generalization of the Moran's formula from Appendix 
IbI It implies that if the IFS /i, . . . ,/jv satisfies ||/j|| < S, then the Hausdorff 
dimension of the attractor admits the following estimate: 'Dy < log Nft/ log 4 , 
where $ is the maximal multiplicity of the continuity partitions for /j. 

Now we claim that as E c — > we have: 5 = max||/j|| < (E c ) rT for some 
a > 0. To see this let us estimate the maximal duration of the avalanche 
f m = vi3x.(i lX )eAxM T (ii x )- This quantity tends to oo as E c — » 0, but not as 
fast as r m (E c , e, A) <~ C\/E c (see H1.2| l. Namely we state that f m ~ C% \og\/E c . 
Actually, if we drop energy 1 to an arbitrary site from a configuration in M, then 
in a finite -E c -independent time all the sites become overcritical. They remain 
overcritical, while the system does not loose a substantial amount of energy. 
During this process the total energy is dissipating in geometric progression with 
an average contraction rate 1 — < 1. So the duration of this stage has 
asymptotic C3 log(l/_E c ). The remaining time to finish the avalanche has a 
smaller asymptotic. 

By m.2\ n(E^. e, A) < Co for small E c . Thus the proof of Theorem implies 
that for certain ^-independent constant C n — kCo for each sequence of C n 
steps in the avalanche process the product of the corresponding S'-matrices will 
have norm < c(e) < 1, which is a uniform estimate in E c -C 1. The number 
of steps in one avalanche grows as logl/J5 c . Therefore 6 < C 4 c(e) c ™ '"s 1 ^ < 
exp(— oTog l/E c ), where a — \C~ X log as was claimed in the estimate. 

Next we claim that 1} < <p(E c ) with (f — o((l/E c ) v ) Vu > 0. In fact, we 
described above the avalanche process for small energies. The first stage is 
finite and contributes only a bounded number of singularity hyperplanes. In its 
seconds stage all of the sites are excited, so there the corresponding number of 
singularity hyperplanes equals the duration. The last stage is shorter of time 
ip = o(logl/E c ), but the number of singularity hyperplanes grows faster, but 
still is bounded by e ^ losN - = o{(l/E c ) v ) (see jOJ for any v > 0. 

Finally D y (e, E c ) < log Nd/ log \ < §^f c - as E c - 0. □ 

Corollary 3. For big £ c >l/e the system [A, F) is not topologically transitive. 

Proof. Suppose (t, x) is a point with the dense orbit in A. We know that for big 
E c the Lebesgue measure of the spatial part y of the attractor A is a positive 
number u> > 0. It follows from Theorem [S] that under iterations with fixed 
excitation sequence t the volume of the spatial part M decreases in geometric 
progression with the number of avalanches. Thus after a finite number of steps 
it becomes less than w. This iteration will be still a finite number of polyhedra, 
so that its closure does not coincide with y. Since it contains all the points 
TT s {F n (t 1 x)) 1 we obtain contradiction. □ 

In the case N = 2 and e > 1/2, the value of E c starting from which 
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%)y(e,E c ) = 2 can be calculated precisely because even one shift of the sloped 
strip mentioned in the above proof overlaps with itself and is sufficient for ob- 
taining an open domain in the attractor. This condition e > 1/2 together with 
E c ^> 1 from the theorem ideologically coincide with the sufficient conditions 
for invertibility of the differentials of avalanche maps (Proposition EJ) . This 
makes an indication of a relation between this invertibility and fractality of the 
attractor in the spirit of Ledrappier- Young formula |LY| . This latter is however 
unappropriate in our situation. 

Note on the usage of the Ledrappier- Young formula. This formula, 
essentially used in |BCK| in the study of the Zhang model, cannot be used for 
the map F : T> — ► T> since this map is never invertible (in loc. cit. it was applied 
to F^ 1 ). In addition to invertibility the Ledrappier- Young theorem is based on 
the SRB-property. For the map F^ 1 this property is equivalent to absolute 
continuity of the stable foliation for F w.r.t. the measure /i. If the measure 
has fractal support this cannot happen. Therefore all formulas based on this 
property may turn to be wrong. We demonstrate this in Example A of 

On the other hand, as we have just shown in the theorem, in thermody- 
namic limit E c — > oo the fractality is lost and so the absolute continuity prop- 
erty is restored (but only for the geometric attractor, the support of SRB- 
measure is smaller!). This however does not help with non-invertibility of the 
factor (E]y, cr~^). Even if we change this factor to invertible two-sided sequences 
(Sat, (Tjv), the system remains non-invertible since not all points of the attractor 
(which can be quite fat) admit negative iterations (Remark 0J. In addition, a 
new negative Lyapunov exponent — log N in the first factor appears and the 
formulas exploited in '!'>( '\\. become completely inadequate. 



6. Examples 

In the examples below we consider the one-dimensional Zhang model with two 
sites, N = 2. 

Example A: A computation shows that for E c > (1 + e)/(l — e) we have six 
domains of continuity [i] x My, i = 1,2. The domains Mxj are given by 



M 11 = {x^[Q,E c ] 2 \x l + l<E c } 

M 12 = {x£ [0, Ec} 2 I xx + 1 > E c and (1 - e){xx + l)/2 + x 2 < E c } (11) 
Mis = {x e [0, E c } 2 | Xl + 1 > E c and (1 - e)(xx + l)/2 + x 2 > E c } 

and the domains M 2 j are symmetric to these. The maps Fij are of the form 
Fij(x) = Lij(x + ei), where 

(l±e\ 2 !-<; " 

y 2 > 2 



Lxi = 



1 
1 



Ll2 — 



e 

l-e 
2 
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Figure 6: The figure shows the physical attractor y = {a,b,c} and the maps 
Fi\y,F2\y for E c = 7/2 and e = 1/2. The arrows on the picture to the left shows 
how the points of y are mapped under Fi , and the picture on the right shows how the 
points are mapped under Fi. 



1 
1 



L12 



l-e 
2 

e 



l-e /l+e\2 
2 V 2 V 



The maps and i^i correspond to avalanches of size 0, the maps F\i and 
i<22 correspond to avalanches of size 1, and the maps F\z and F23 correspond 
to avalanches of size 2. 



It was discovered in jBCKj that the physical attractor y has the following 
simple structure: 



y 



1 



l + e N ,1 + e 1 



l-e 2 



e' 1 



1 + e i+_f O 
1 - e' 1 - e'J 



for S [frf , jz^} We denote these points by a, b, c so that y = {a, b, c}. The 
maps -Fi|y and i^ly are permutations of 3^ 



6 c 
a 



Figure El shows the physical attractor y and the maps -Fi|y and i^y for 
£ c = 7/2 and e = 1/2. We construct a partition 7£ = {Rx, R 6 } of E J x F 

by 

Rx = [1] x {a} R 2 = [1] x {b} R 3 = [1] x {c} 
i? 4 = [2] x {a} R 5 = [2] x {&} i? 6 = [2] x {c} ' 

We let ^4 = || ay || be the 6x6 matrix where ay = 1 if F(Ri) n 7^ and 
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dij = otherwise. It is easy to verify that 

1 0" 
1 

1 



1 

It is clear that gn : A — > is a topological conjugancy of the maps F\a 
and a\. The matrix A is transitive and Sp(^4) = { — 1,-1,0,0,0,2}. Hence 
^topC^U) =log2. 

If \i is the SRB- measure on A, with (tt u )^ij. being the uniform Bernoulli 
measure on S^, then (giz)*n is the Perry measure on With respect to this 
measure it is easy to see that the average avalanche size is so = 1. It then follows 
from that the sum of the negative Lyapunov-exponents is xi + X2 — l°g e - F° r 
instance we see that for E c = 11/2 and e = 2/3 we have h^F]^) > \xi\ + | X2 1 - 
So even though Fi\y and F 2 \y are invertible, the Ruelle inequality (and therefore 
the Pesin formula) cannot be reversed (though this was argued in p3CK ). 

We also remark that for E c = (1 + e)/(l — e) we have y C S(F), so the 
standard construction of SRB-measure will fail in this case. Still there is clearly 
a natural invariant measure. 

The example with a trivial physical attractor can be generalized to all N for 
d = 1. Then y consists of N + 1 points Zo, Zi, . . . , Zn given by 

__ 1 + e 
1 — e 

where eo = and ei, . . . , ejv is the standard basis in K . 

Example B: For E c = e = 1/3 the dynamics is very simple. The map Fi 
has two domains of continuity: Mu — {(x,y) £ M\y < — 2x/3 + 1/3} and 
M12 = M \ Mu. The domains of continuity for F2 are given by symmetry. The 
maps are given by ^11(2;) = Lu(x + ei) and Fi2(x) = Lyi(x + ei), where 

2 3" 
2 3 ' 

We see that all four Fij are mappings to the diagonal line in M. In fact, 
Fn(Mu) = F 2 i(M 2 i) = [ Pl ,p 2 ], where Pl = (2/9,2/9) and p 2 = (1/3,1/3). 
The interval is contained in M 12 PI M 22 since the two lines of singularity intersect 
the diagonal in the point (1/5, 1/5). The images of A/ 12 and M 22 also coincide 
and is an interval [^1,733] C [^1,^2], where p% = (22/81,22/81). This shows 
that singularities are removable after one iteration. It is easy to see that for all 
t £ the dynamics will contract to the fixed point P = (4/17,4/17). Hence 
the attractor of the system is A = x {P}. The dynamics on the attractor 
is of course conjugated to <r^. 
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L u = - 
9 



2 3 
2 3 



and L 12 = — 
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The example can easily be extended to a neighborhood of (1/3, 1/3) in the 
(e, -E c )-plane. In the region 

r 1 - 2e + 13e 2 e(7e 2 - 6e + 3) 1 ^ _ ^ . r e 1 - 2e + 5e 2 1 
maX U + 12 e -13^' 763-10^ + 7,-4/ " ^ " mm tl=? 1 + 4 e - 5e 2 I 

we have the same avalanches as for E c = e = 1/3. This region is shown in 
Figure For all points in the region the maps depend continuously on e and 
the lines of singularity depend continuously on e and E c . The condition for an 
atomic spacial attractor is that the images of the domains do not intersect the 
singularities. For E c = e = 1/3 the images are bounded away from the lines 
of singularity, and hence there exists an open neighborhood of (1/3, 1/3) in the 
(e, .E c )-plane where the same holds, i.e. the attractor is of the form x {P} 
(where P is a point in M) and the dynamics is conjugated to cr^. 

Example C: Let us consider E c — 1/3 and e = 1/2. In this case there are 
28 domains of continuity and 28 corresponding maps. A computer program 
is written to compute the sets U n . The program uses exact calculations of 
the edges of the polygons that make up U n , and hence it can be used to give 
rigorous "proof by computer" of removability of singularities. By using the 
program we obtain that singularities are removable. In fact C/5 l~l S = 0. The 
set U$ consists of 13 connected components. Figure |S] shows the set U§ and 
the lines of singularity, and Figure is a schematic illustration of how these 
connected components are situated with respect to the lines of singularity. 
The intersection of the connected components of U5 with y are denoted by 
Yi, . . . , Y13. Then we construct the partition 1Z = {[i]xYj}, and enumerate the 
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Figure 8: Shows the set U 5 for N = 2, E c = 1/3 and e = 1/2. The points on the 
attractor are magnified in order to make them visible in the figure, and hence it looks 
as if they intersect singularities, but in fact they do not. 




Figure 9: Shows how the 13 spatial partition elements are situated with respect to 
the lines of singularity. 
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Figure 10: Shows the matrix A for the coding of Example C. Black squares are 1-s 
and white squares are 0-s. 



elements so that 1Z — {R\, ■ ■ ■ , R 2 6}, where 

R 1 = [l]xY 1 R 3 = [l]xY 2 ... R 2 5 = [l}xY 13 
R 2 = [2]xY 1 R i = [2]xY 2 ... R 26 = [2] x Y 13 ' 

We construct the 26 x 26 matrix A = ||ciy || by letting = 1 if F{Ri) n Rj 7^ 0, 
and a,ij = otherwise. After making the computations, the matrix A becomes 
as shown in Figure 1101 The black squares represent ones and white squares 
represent zeros. Direct computation shows that the matrix A is transitive. 

Since singularities are removable the map g-jz : A — > is an avalanche con- 
jugancy between and a\. The SRB-measure projects to the Perry measure 
on £j[, so it is possible to calculate properties such as average avalanche size. 
In this example a computation gives sq = 123/17. The spectral radius of the 
matrix A is 2, and hence th op (<j\) = log 2. 

Example D: In the previous examples singularities are removable. This is 
however not always the case. FigurelTTlshows the set U 2 q for N — 2, E c — 7 and 
e = 1/2. This is an example where singularities are non-removable. We will in 
the following show that singularities are non-removable for E c = (3 + e)/(l — e). 

Since (3 + e)/(l — e) > (l + e)/(l — e), the domains of continuity are given by 
the formulas presented in Example A. Take the point p = (E c , E c — 1) S M\ 3 . 
See FigureEl Observe that p lies on the horizontal line x\ = E c — 1 and hence 
p G S(F). Clearly p is in the interior of M13, so 

w-MBS-'l-K-^-f'-M' 

Denote q% :— F\{p) and observe the it lies on the singularity line x 2 = E c — 1. 
On the other hand qi is in the interior of M 2 \ 1 so ^2(91) = p — (1, 1). Denote 
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Figure 11: Shows the set U20 for E c — 7 and e = 1/2. In this example singularities 
are non-removable. 




Figure 12: Illustration of the fact that singularities are non-removable for E c = 7 and 
e = 1/2. 
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1i '■= ^2(91)- This point also lies in the interior of M 2 i. Let q 3 := F 2 {q2) = 
p — (1,0). It is clear that Fn(q 3 ) = p, and in this sense p is a periodic point. 
However, i*i(g 3 ) is not well defined since q 3 £ 9M 12 n 9M 13 . 

In the following we let (z\, . . . ,z n ) denote the open convex polygon with 
edges Zi, . . . , z n . Define Bg(p) — (p 7 p + (0, S),p+ (—5, S),p+ (—5, 0)). For small 
5 > we have Bg(p) £ M13, and hence 

Fi(B s (p)) = F 13 {B s (p)) = (qi,q x + 5a, qi + 5b, q x + 5c) , 

where 



1-e \ , /l-4e-e 2 lit \ / A + e, 2 1-e 



b = .«= -(^) - 



2 y ' 2 



-6 



The polygon Fi(Bs(p)) intersects the singularity line x\ = E c — 1. See Figure 
1121 A simple computation shows that 

Fx{B S (p))r\Mu = (qi,qi+5a,qi+5b,qi+5c) , 

where 

(o.(^)0- 

It then follows from the above discussion that 

Fx o FiiF^Bsip)) n Af u ) = <P,P + 5a 1 ,p + 56, p + (5c) . 

It is then easy to verify that for all 5 > there is 7 > such that B 1 {p) C 
tt s o F i {Y,+ x B§(p)), So for each n £ N there is 5„ > with Bj n (p) C [/„. The 
image of B$ n (p) under F\ intersects singularities, and its closure contains the 
point qi, which thus is an essential singularity. Clearly the points p, q 2 and q 3 
are also essential singularities. 



7. Statistical properties 

In order to evaluate the entropy and Lyapunov spectrum of the physical model in 
the thermodynamic limit we need to derive several estimates for the asymptotic 
behavior of observables like avalanche size, avalanche duration and "waiting- 
time" between avalanches. The results are derived using only the uniform 
Bernoulli measure on S^, and hence hold for any SRB-measure and for time- 
averages. 

In ^17.31 we define the thermodynamic limit as the double limit E c — > 00, 
L = \/~N — > 00, contrary to jBCKj . where the thermodynamic limit is defined 
as the limit L — ► 00 only. This is important as the quasi-classical limit since 
equivalently means a fixed energy, but the energy quantum of Section \T\5 — > 0. 
As E c — > 00 we must make a scaling of time in the physical model. Otherwise 
the influx of energy to the system will go to zero. With this new scaling we 
show that the entropy goes to zero and the Lyapunov spectrum is collapsing. 

We do not provide strict mathematical proofs, but still think important to 
include the discussion of our results from the physical point of view. 
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7.1. Statistics of observables 

Let t be the coordinate measuring the duration of avalanche and let u corre- 
spond to the interval between avalanches (minimal value 1). We will also study 
the observable s - the avalanche size (defined in HI .3(1 . While in the first case 
we consider only actual avalanches, so that r > 0, in the second we make dis- 
tinction between sq - all avalanches including the trivial case of under-critical 
state (s = 0) and s + - the actual avalanches, so that s + > 0. 

The reason for introducing two different avalanche size observables is the 
following: sq plays a crucial role in mathematical investigation of the model 
(see H7.4[l . while s+ is important from physical perspective. In we will see 
that physical observables should allow a thermodynamic limit. 

Denote by f, w, so, s+ the corresponding mean time-average quantities, each 
of which is a function on the space- factor M and is defined as follows: 

a(x) = / lim - Vcr(/ l (t,a;)) d/x Bor 

/v+ fc-+oo K L — ' 

with fj, Ba[ being the Bernoulli measure on the one-sided shifts (by Birkhoff cr- 
godic theorem the time-average limit exists almost everywhere and is measur- 
able). This function is invariant in the sense: a(x) = jj Y^i=i &(f{h x ))- 

Whenever the system is ergodic with respect to an invariant measure /i the 
function d is constant /x-a.e. and equals the space (ensemble) average 




If the system has a unique invariant SRB-measure the function a is constant 
/i Bcr x /i Lob -a.e. and equals the space (ensemble) average 

(cr) = J a cfyi SRB . 

S+xM 

We will need the maximal values of these observables in a sequel, which 
we denote by r max , cj max and s max = maxso — maxs+ respectively. We shall 
calculate their asymptotics in L and E c . 

Denote by tp\ ~ ip2 the asymptotic equivalence relation meaning that the 
ratio ifi/if2 has subexponential grow/decay. We denote the equivalence by 
when the limit of the ratio is 1. 

Lemma 21. For E c 3> 1 the maximal avalanche time and size have the asymp- 
totics: T max ~ L 7t , u max = L d E c , s max - L d+ls , where j T = j g = 1 for d = 1 
and 1 < 7 r , 7 S < d for d > 1 . 

Proof. Consider at first the simple case d = 1. The maximum avalanche 
duration and size are achieved when all sites contribute to the avalanche, i.e. 
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their energies are sufficiently big and there is one site with energy greater than 
E c — 1 to initiate the avalanche. Actually, if some site has small energy, it will 
serve as a boundary and the avalanche wave reflects from it (to be explained 
below). 

Assuming E c > we know from Lemma that in the avalanche process 
each site, whenever overcritical, relaxes until in the next step it receives a suf- 
ficient portion of energy to become overcritical and relax etc. In other words, 
the sites blink, being under- and over-critical in turn. But in this process they 
make overcritical their neighbors and the process propagates as a wave, with 
only difference that its front excites new sites, while in the traversed region there 
remain blinking overcritical sites. 

This wave spreads along the interval B\ =[l,L]cZ towards its boundary 
and then it reflects from it, bearing now relaxation. In fact, as the front wave 
reaches the boundary it losses a substantial part of the energy on the boundary 
sites, which thus cannot be recovered and remain undercritical for the rest of 
this avalanche. They influence their neighbors to stop being critical and so 
forth. Thus we obtain the reflected wave that, in contrast with the first one, 
turns overcritical sites into relaxed. When the wave hits itself, the avalanche 
process stops. 

It is clear that the duration of this avalanche is r max sa L. The number of 
involved sites corresponds to the area of the triangle with a side B\ — [1,L] 
and height L/2 (recall that each site is overcritical only half time of its blinking 
period), whence s max « i 2 /4. 

Consider now the case of dimension d > 1. Here the scenario of maximal 
avalanche is more complicated and consists of three stages (the proof is similar 
to the case d = 1, but quite lengthy and will be suppressed). Again the maxi- 
mum avalanche duration and size are achieved when all sites contribute to the 
avalanche, though now if some isolated site has smaller energy it serves as a 
boundary only once but then on the next several waves it receives the required 
portion of energy and follows the general scheme of motion. 

At the first stage a site is excited and it initiates the rhombus-shape wave 
(a cube in the Manhattan metric, see Figure IT3)) that spreads to the boundary 
of the cube Bf — [l,L] d C Z d (in time ~ L/2 if the center of the rhombus 
is placed near the center of the cube). The second stage begins as the wave 
reaches the boundary face and reflects from it (See Figure H^ . The reflected 
wave is almost momentary overthrown by the coming overcritical wave, which 
again reflects from the boundary, come now deeper into the interior of the cube, 
but is overthrown too etc. If one looks along the boundary face, the reflected 
wave travel along towards a vertex with preserved form (like a soliton) and then 
disappears into this vertex (there occur strong interactions with other waves in 
this corner). But if one looks into the perpendicular direction, the collection of 
reflected waves oscillates (each reflected wave enters deeper and deeper into the 
cube) contributing to the avalanche duration the sum 1 + 2 + 3 + . . . , which stops 
with the end of the second stage (we do not specify the sum precisely because 
after some oscillations the wave front becomes more and more eroded by the 
interactions between overcritical and relaxing waves; this impairs the sum and 
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Figure 13: The picture is a "snapshot" of the lattice A for d = 2, L — 10, E c = 7 
and e = 1/2. A single site in the marginally stable configuration has been exited, 
and a great avalanche is unfolding in a rhombus shape. This is the first stage of this 
avalanche. The white squares are overcritical sites, the gray squares are sites with 
energy just bellow E c and the black squares have energy approximately equal eE c . 



decreases the exponent, but not too drastically). 

The third stage begins as the main body of the overcritical sites becomes 
disconnected and the avalanches behaves like worms crawling along the high 
energy fractal-like collection of states. We illustrate this in Figure ED From 
the description of the maximal avalanche process it is clear that the asymptotic 
exponents 7 T ,7s do not depend on the energy E c ^> 1. Let us denote (the 
sequences are increasing, so the limits exist): 

,. logT max logs max 
7 T = Inn — ■ and a + j s — lim 



L->oo logL L^oo logL 

Duration of the first stage of avalanche is ~ L. The above arguments show 
that the exponent j T of the second stage is > 1 . The last stage can only increase 
it: 7r > It- To see that j T < d we note that in average the number of critical 
sites on the boundary is about kl ~ L 7 " with < 7 K < d. Thus we have a 
constant flow of energy out of the system with the average speed > ^j-klE c , 
and the inequality E tot < L d E c proves the claim. 

To estimate the maximal avalanche size exponent 7 S consider again the sec- 
ond stage of the above scenario. The number of involved sites corresponds to 
the volume of the prism IIj, over the cube with height £ w ^T max ~ L 1 ^ , i.e. 

L d +< ~ Vol d+1 (n L ) = -J— ■ Vo\ d (B d L ) ~ L d+ <. 

Thus 7 S > 7s = j' T > 1. Inequality 7 S < d follows from the inequality from 
above for the duration exponent 7 T . 

The maximal value of the waiting time is obvious. □ 
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Figure 14: The picture shows a "snapshot" of the second stage of the avalanche 
shown in Figure ED We can see the well-shaped (soliton-like) waves of energy near 
the boundary and observe how their form begins being eroded near the vertices. 




Figure 15: The picture shows a "snapshot" of the third stage of the avalanche shown 
in Figure 1131 The energy configuration has a fractal structure where the avalanche 
can sustain in regions of high energy. 
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7.2. Asymptotic of the statistical data 

Now we can study statistics of the avalanche data asymptotically (as for the 
thermodynamic limit). 

We will need the following technical statement (informal only for E c > 1): 

Lemma 22. Almost every (w.r.t. a random excitation sequence) spatial trajec- 
tory returns to the cube Bq = (E c — 1, E C ] N . 

Proof. It follows from Proposition 0] that K = {x G M \ 3i : Xi S (E c - 1; E c ]} 
is a return set, i.e. every trajectory F n (t,x) meets x K. Partition the 
spatial part y of the attractor according to the hyperplanes collection H = 
^>i,men{xi = E c -m}. 

Each such a part can be shifted by an excitation sequence to the cube Bq. 
The probability of all such sequences (where the avalanche starts from the set 
Bq of maximal energy in all sites) is positive. Let p > denote the minimum of 
these probabilities over the finite set of all partition elements of Y by 7i. Then 
the probability of not entering Bo in k successive avalanches is less than (1— p) k . 

Therefore since the number of avalanches tend to infinity as we iterate the 
dynamics, the measure of trajectories staying away from Bq is zero. □ 

Theorem 23. We have the following asymptotic estimates valid as E c — > oo 
(and N 3> 1 fixed) or L = \/~N — > oo (and E c 3> 1 fixed): 

(i) f ~£^; 

(ii) Q ~ E c ; 

(hi) Sq ~ L d +^/E c ; 
(iv) s+ - L d +^ , 

where "y T , 7 S are the same exponents as in Lemma Wl\ (thus 7^ = lim }° s< ^ — 0). 

Proof. The maximal avalanche size is achieved for a certain configuration of 
states V C M, which we can bound as follows: Ui C V C Ud, where 

Uj = {x £ M I Vi : (1 - j^)E c < x t < E c and 3i : x lQ > E c - 1}. 

Denote also Uj = {x e M \ Mi : (l - j 1 ^i)E c < x l < E c } D Uj. 

To estimate the measure of the sites leading to the maximal avalanche we 
consider preimages of x Uj under the map F. It is clear that one needs 
k e G [jzr, different backwards iterations F~ 1b , s — 1, . . . , k e , to cover the 
spatial attractor y. Since the measure p is F-invariant, we get for its 7r s -push- 
forward: p s {Uj) « Pj/k e , which is ^-independent. Thus we get the same 
exponent (d + j s ) for s+ as for s max - 

The same arguments yield the asymptotic of f . 

To obtain the asymptotic of w in L we note that since the amount of lost 
energy is < CL 1k+s E c (where j K < d is the quantity from the proof of Lemma 
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Figure 16: The figure illustrates how the continuity domains are mapped under i*\ 
and F 2 for N = 2 and E c > e/(l - e). 




Figure 17: The figure illustrates how the faces of the continuity domains are mapped 
under Fi, F2 and F3 for N = 3 and E c > e/(l— e). The point of view is at (E c , E c , E c ). 



I21|) . the average remained energy in a site of configuration obtained from a 
maximal one after an avalanche is (L d E c — C , L 7re+<5 E c ) / ' L d ss E c . Thus the 
waiting time does not grow with L. 

The asymptotic of Q in E c is quite different: If N is fixed but E c grows, 
then any state from dM becomes at distance On ■ E c after some relatively small 
number of iterations. 

Let us first demonstrate the idea in the simple case N = 2. For critical 
energy E c > e/(l — e) the picture of avalanches is shown on Figure El We see 
that in a few steps of the dynamics the configuration becomes far from dM, 
i.e. it strongly contracts in all directions. Actually, it is possible to imagine the 
situations when the point is mapped to the vertical strip and then is shifted 
horizontally for a long time by excitations of the first site, but probability of 
this event exponentially goes to zero as E c — > 00. Thus in a relatively short time 
the point from dM is mapped into the square [0, 92E C ] 2 , where the constant 82 
is _E c -independent. To achieve the boundary dM again it needs ~ (1 — 6>2)2_E C 
random excitations. 

The similar picture happens for N — 3, see Figure El in the general case 
Theorem insures that after some (few) number of steps we get strong con- 
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traction, so that the point becomes in the cube [0,9nE c ) n , i.e. far from the 
boundary dM. Thus we need ~ (1 — 9^)NE C excitations to make it overcritical 
and this implies the claim. Note that the asymptotic for Q is not for the double 
limit E c — > oo, L — > oo, but for two partial limits only. 

To obtain the estimate for so we need to estimate the conditional measure 
Hs(Uj\Uj), which coincides with the probability that a randomly chosen config- 
uration x G Uj and site i satisfy: Xi > E c — 1. This probability is ~ b±/E c . 
Thus ju(S+ x V) ~ a e /E c and the mean avalanche size is (sq) ~ L d+ls / E c . 

This is however the space-average (the arguments below work also for s + , 
f). We would obtain the same for the time-average of Lebesgue a. e. -initial con- 
dition if we have an SRB-measure, or for fi-a.e. if we have an ergodic measure 
/z. But we cannot guarantee existence of an SRB-measure or an ergodic mea- 
sure. However, for a non-ergodic measure, we can decompose A into ergodic 
components, where the Birkhoff theorem works. By Lemma 1221 each ergodic 
component (with non-trivial contribution) intersects in its spatial part the top- 
energy cube B and so for each of it the same asymptotic of space-averages 
holds with universal exponents but maybe different coefficients. Therefore we 
obtain the required asymptotic for sq. 

Another way to get the last asymptotic is via the formula sq — - '-^tj^ ■ D 

For d — 1, E c £ [jzri we already have f ~ N = L 1 as the theorem 
states, but so ~ N = L 1 , N — > oo, which shows in the last respect the critical 
energy E c — 2/(1 — e) is small. 

In |BCK| the estimate 7 T > 1 was predicted for all d, while this is a feature 
of the cases d > 1. In the latter cases the analytic calculation of exact values of 
exponents 7 r ,7s is a difficult problem. 

Remark 9. The difference between cases d — 1 and d > 1 demonstrated in 
the theorem is known in the physical literature. The former case is usually 
considered as the trivial SOC-model. 

7.3. Thermodynamic limit 

By thermodynamic limit of an observable (f> we understand the double limit 

[4>]oo = lim (j) 

L — >oo 
E c — *oo 

if it exists. It is assumed that for physical observables this limit exists. In |BCK| 
only limit L — > oo was considered, though then the value of energy E c could 
serve as an essential parameter, which is not desirable in the SOC-paradigm. 
However it was suggested there that consideration of E c — > oo can be helpful. 

As an example of non-physical observable we expose sq (in i j7.1l we called it 
mathematically relevant): The double limit does not exists because the repeated 
limits are different: 

= lim lim s ^ lim lim s = +oo. 

L — >oo E c — >oo E c — >oo L — >oo 
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But s+ and r are good physical observables, for the thermodynamic limits exist: 

[f]oo = oo, [s+]oo = oo. 
From fuse iteration arguments in the topological case) and H4.2l we obtain: 
[^to P (/)]oo = [htopiF)/^^ = and [hto P (F)]oo = oo. 

But with lo the situation is different because the proof (rather than the vague 
statement of part (ii)) of Theorem 1231 implies: lim^ ^ lim^^oo tD — oo, while 
Ym\E c ^oo hniL-+oo w is finite. Since to is definitely physically relevant observable 
one needs the following reparametrization: to t— > lo/E Ci which corresponds to 
contraction of the waiting time via the following ansatz: 

We let the energy quantum added at a unit time to the system equal 5 = h 
(instead of 1 as before), but speed up time in the waiting intervals respectively: 
E c = Eq/H, uj ncw = ujfi. Then the thermodynamic limit (the space part M can 
be quantized similarly via L — [l/h] with I a finite length) corresponds to the 
quasi-classical limit h — ► 0. 

The duration of avalanche was suppressed in our definition of dynamics to 
length one. So if we want to find the entropy of the physical system, where 
each step of avalanche has time-duration one, we should multiply it by the 
probability of dropping energy into the system. For every trajectory this equals 
^Xf. This ratio behaves differently as L — > oo or E c — > oo, so we need the 
reparametrization described above. 

In M4.2I we calculated the entropy of the "return" Zhang model h t o P (F) = 
log N (this was proved almost surely, but even with the possible entropy growth 
for some exceptional parameters the thermodynamic limit below is unaltered). 
But after reparametrization it changes. Denoting by h zhans the entropy of the 
reparametrized system we get: 



•'new \ 

_1_ ,= / max * 



hf^ = VK) • <^^>, C anS = d-\o g L.{- 
This implies: 

[^ hang ]oo = [Cp ang ]oo = 0. 

Therefore the expanding property is lost in the thermodynamic limit for the 
original physical system, as was already noticed in BCK for a bit different 
situation. 

Remark 10. Notice that in the reparametrized system uj new <C f, which is 
counter-intuitive for certain SOC-examples (sandpile, earthquakes etc, where 
one expects u> 3> f). This indicates that the Zhang system should be modified 
by introducing the local contraction of time depending on the avalanche size or 
speed. We will not consider such gradient-type models here. 
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7.4. Lyapunov spectrum 

In M.HI we showed that the Zhang model is hyperbolic with one positive exponent 
and the remainder of the spectrum negative. This hyperbolicity is lost in the 
thermodynamic limit. 

Proposition 24. For E c > e/(l — e) we have 2j=i Xi = *o loge. 

Proof. We cannot ensure the existence of an SRB-measure, so both the Lya- 
punov exponents and so should be seen as functions on M . From the general 
theory of Lyapunov exponents we know that 

N „ 71-1 



EX7(%)= lim - VlogdetTfT t x) . 
n— >oc n — ' 



n — >oo Ji 
i=l t=l 



From Proposition we know that the formula det(Ly) = e 3ij holds for E c > 
e/(l — e). Hence 



N N 



0) = ^2 + Xi {t,x)d(i Bet 

i=\ i=l Jt, n 

f N 

J Tit 



r i n_1 

/ lim - Vlog(detT)(F*(t,a;))d/x E 
;„ + n^oo n £ — ' 

N t=0 
f 1 

/ lim — y s(t, x) dfj, B<!I loge 
n^oo n L — ' 



t=0 

s (a;) loge. 



Corollary 4. |x | = 77 log \ — ► as _E C — > oo } om£ |x | — > 00 as L — > 00. 



□ 



Thus we should study differently the following cases: 

1. E c — ► 00, but L (and d) fixed. Since UruE c ->oo So = 0, the negative part 
of the Lyapunov spectrum collapses: lim^; < .^ 00 Xi — f°r every 1 < i < N. 
The hyperbolicity is lost, but the positive exponent Xo = log TV survives. In 
particular, the entropy does not collapses. 

2. L — ► 00, but i£ c ^> 1 fixed. Here only a bounded piece Xii---iXk °f 
the Lyapunov spectrum collapses, k — const. But the number of elements of 
this spectrum grows and in average \x~\ — * 00. In particular, |x~|max — » 00 as 
L — ► 00. Again, the positive exponent xo and entropy are preserved and though 
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we loose hyperbolicity there are many non-degenerate Oscelledec modes. More- 
over they prevail over collapsing modes and so essentially the hyperbolicity is 
preserved as well. 

3. Reparametrized model. This was introduced in H7.3I and require the 
renormalization: multiplication of waiting time by the function along the 
trajectory. In this case the Lyapunov spectrum collapses to zero in any limit 
E c — > oo and L — ► oo and the hyperbolicity is completely lost. 

Thus the exponential grow of the statistics is suppressed and we can observe 
power law statistic as is the basic idea of SOC-phenomenon. The corresponding 
SOC-exponents are related to the asymptotic of the Lyapunov spectrum (as 
discussed in ffl ( ' K ) . but are difficult to calculate analytically. 

8. Conclusion 

It is of importance to the paradigm of SOC to understand the mechanisms 
behind the behavior one observes numerically in the Zhang model and other 
sandpile models, and the aim of this paper is to provide a first step to a rigorous 
mathematical understanding of the dynamics of Zhang model. 

Due to the singularities and non-invertibility of the model there existed very 
few applicable results, and hence we had to modify known results and develop 
some new methods in order to describe the dynamical properties of the model. 

The result of this work is that the singularities play a modest role in the sense 
they do not change the main dynamical characteristics. However the effects of 
the singularities can be seen in the rich fractal structure of the spacial attractor. 

Our analysis allows to take the thermodynamic limit of the main dynami- 
cal quantities, showing that the entropy vanishes and the Lyapunov-spectrum 
collapses after re-scaling the model. From a physical point of view this is in- 
teresting. Typically chaotic dynamics (positive entropy and positive Lyapunov 
exponents) is an indication of exponential speeds of mixing (short decay of cor- 
relations) |Ba| . which is not compatible with the power-law statistics of the 
SOC-hypothcsis. The loss of hyperbolicity in the thermodynamic limit hence 
supports the critical behavior observed numerically Q], |7], |GD| . 

To conclude: We have shown that the Zhang model is a chaotic hyperbolic 
dynamical system, where all the entropy is produced by the random driving of 
the system and, due to singularities, the orbit structure is richer than for a topo- 
logical Markov chain. The hyperbolicity indicates that under weak conditions 
there exists an SRB- measure (self-organization). In the thermodynamic limit 
the hyperbolicity is lost and we may expect power-law statistics (criticality) . In 
practice the systems that are studied have finite size and finite critical energy. 
Hence they are chaotic, but with small entropy. The SOC-hypothesis is that 
these weakly chaotic systems have SRB-measures with the rates of convergence 
to these measures being exponential but slow compared to a unit step in an 
avalanche, causing the prevailing of power-law statistics. 
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Finally note that modifications of the Zhang model are possible. For in- 
stance, different amounts of energy Si can be added to different sites i G A in 
the excitation process. This corresponds to the rectangular form of A C Z d 
and uniform quantum of energy 5 = H. We can also consider other spacial con- 
figurations A. In this way nearly Zhang models considered in Section 14.21 are 
natural. Another approach is to use a random amount of energy, i.e. stochastic 
S : A — > R + (this idea was used in the plasma physics |KK| L which can also 
be presented as a skew product with piecewise affine fibers, but now over a 
solenoidal system. In all these theories the ideas from the present paper work 
well (though the thermodynamic limit will be sensitive to the form of A). 

A Topological entropy of piecewise affine maps 

A.l. The Buzzi theorem and its generalization 

The statement as it is done in |B1| does not apply to our situation. In B2 Buzzi 
noted that it extends to isometries and contractions. In fact, the assertion holds 
always, but since we cannot make a simple reference we write an adapted proof 
for the convenience of the reader. 

Theorem 25. Let X C M. d be a bounded polytope and f : X — > X a piecewise 
affine map. Then 

H sins (f)<X+(f)+H mun (f), 

where 

X + (f)= lim sup - max log II A^/™ || . 

Proof. Let V = {Pi} be the continuity partition and /j := f\p i be affine maps. 
Fix e > and T = T(e) > f log(Vd + 1) such that for n > T we have: 

mult(7>") < exp((ff mu]t (/) + e)n), \\Adf n \\ < exp((A+(/) + e)n). 

Take r = r(e) to be compatible with the partition V T (for f T ), i.e. any r-ball 
intersects maximally mult('P T ) partition elements. 

We will prove that each non-empty cylinder C(a) — [P ao . . . P aiT1 ] of length 
|a| = IT can be partitioned into a collection Q(a) = {W} satisfying the following 
properties: 

1- E|a|=iT card (Q( a )) <C exp((A+(/) + ff mult (/)+3e)ZT) 
2. diam(/„ l!r _ 1 o...o/ O0 (W))<r. 

Let us prove this claim by induction assuming it holds for some I > 0. The 
base of induction is obvious and C'o is the minimal cardinality of an r/2-ball 
cover. 

Take a partition element W £ Q(a) that is used to cover the cylinder 
[P ao . . . P aiT1 ]. By the induction hypothesis it has diameter less than r, so 
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it can be continued to cover a non-empty cylinder of length (I + 1)T in at most 
mult(P T ) ways. So to cover the cylinders [P ao . . . P aiT _ 1 Pb Q ■ ■ ■ Pb T -i] wc ma ke 
a division of W: 

7 

W=\J\Vl ,7< multCP T ). 

i=l 

Let 

Wl' = f aTl _ 1 o...of ao (Wl) 

and 

Wl" = f bT _ 1 o...of bo {Wl l ). 

By the assumption diam(W/') < r for alH = 1, . . . , 7, but the sets W(" may have 
greater diameter. We need to divide the sets W-" so that they have diameter 
less than r, and then pull this refinement back to the partition of sets W-. 

Let L be the differential of fb T -i ' ' ' h a on W" . We can assume that 
L is symmetric and take {e^} to be a basis of eigenvectors corresponding to 
eigenvalues Ai, . . . , A^. Let be a basis in the vector subspace corresponding 
to W-". We can choose this basis to be orthonormal and triangular with respect 
to {e k }. Divide W-" by the hyperplanes 

^(x) d = (vj,x) =P-J=, P eZ . j = l,...,d. 

This defines cells W of diameter less than r. Since i>j(W(") = ipj(L(W")) has 
diam < |A,|r, the number of cells W needed to cover W-" is less than or equal 
to 

(Vd+ l) d |A!|+ . . . |A d |+ < (Vd+ l) d \\Mp\\ < cxp ((A+(/) + 2e)T) , 

where |Ai| + = max{|A,|, 1}. 

Therefore the total cardinality of the new partition is less than or equal to 

mult(7> T ) exp ((A+(/) + 2e)T) cxp ((A+(/) + H mult (f) + 3e)lT) 

< exp ((A+(/) + H mult (f) + 3e)(i + 1)T) . 

This proves the statement. □ 

The theorem holds as well for most degenerate piece- wise affine systems, but 
there can be problems with H mu \ t (f). Namely the latter is not defined if the 
image of continuity domain contains a boundary face of a continuity domain. 
But if we assume the image and the faces always meet transversally, no problems 
occur and the above theorem applies literally. 

In degenerate cases of the Zhang model, the above requirement holds for 
most parameters. For instance, if N — 2 and e = 0, then all E c ^ Z satisfy the 
request. For integer E c the theory fails, but one can look just to the whole image 
set, which has dimension < N: its Poincare return map is piece-wise affine and 
it satisfies the requirements. In the above example N — 2, e = the dynamics 
is confined to two one dimensional lines, whence the multiplicity entropy is zero 
(by dimensional reasons) and h top (F) — log 2 in this case. 
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A. 2. Entropy of conformal piecewise affine skew-products 

We say that an affine map is conformal modulo degenerations if the image on 
some subspace transversal to the kernel is mapped conformally to its image. 
A piecewise affine map is said to be conformal modulo degenerations if all its 
affine components are conformal modulo degenerations. We will assume that 
degenerations satisfy the transversality requirement of lA.ll 

In |KR1) we noticed that -ff mu i t = for piece-wise affine conformal maps. 
This easily extends to allow degenerations. Now we consider a more general 
situation of skew-product systems of Zhang's type. 

Theorem 26. Let fi'.X^Xbe piece-wise affine non-strictly contracting and 
conformal modulo degenerations, i = 0, . . . , N— 1 . Define F : E^- xl^ E^, x X 
by the formula F(t, x) = (<7^-(t), ft (x)), t — toti ■ ■ • G E^y, x G X C W l . Then 
we have: htop(F) = logN (so that a-posteriori the variational principle holds). 

Remark 11. For N — 2 and e = the affine components have rank 1, and 
hence Theorem \2b\ shows h top (F) = \ogN. The same holds for E c — e — 1/3. 

Proof. It follows from Theorem 123 that it suffices to prove that H mu \ t (F) = 0. 

To achieve the desired equality note that preimages of the time-like singu- 
larity planes {t = n/N} never intersect under inverse iterations of F(t,x) — 
(Nt mod l,f[ Nt] (x)). Thus H nmit (F) = sup t if mu it(/t)< We will prove that 
flnmit(/t) = for all t e E+ . 

A piecewise affine map can be considered as an ordered triple (X,V,f), 
where X is a polytope in R d , V = {Pi} is a partition of X made up of pairwise 
disjoined polytopes (with certain faces of boundary included, so that the whole 
boundary is distributed between polytopes) and fi := f\p t : Pi — >■ X are affine 
maps. We let X' = U tut(Pi) and Sing(/) = X \ X' . 

For a piecewise affine map (X, V, f) and a point x € X construct a piecewise 
affine map (X X ,V X , f x ), called the differential of / at x, by letting 

1. X x = {y £ R d | 3e > Os.t.Ve € (0, e ) : x + ey e X} C M d is the tangent 
cone to X. 

2. V x is the partition of X x consisting of non-empty sets 

Px = {V G R d | 3e > Os.t.Ve G (0, e ) : x + ey G P} , 
where P G V . 
3- f x 'Px~* Xf( x ) is the collection of maps 

Lun r-,, f: - r -' !l> Hl;l - P xeVx . 

e^0+ e 

Consider iterated differentials and denote f Xl ,...,x n '■= (• • • (fxi)x 2 ■ ■ ■ )x„- 
Note that (f?) x = (/i n _J/»(x) ° • ■ ■ (Mf^x) ° (ft )x- 
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Proposition 27. There exists a constant C £ K+ such that for any subspace 
W C R d and any (xi,... ,x r ) £ X x W r ~ 1 we have: 

mult((/ t "|^) Xl; ... iXr ) < C ( sup sup mult(/ t |y) yi ,...,j /r+1 ) , Vn > 0, 

where the collection of points (3/1, ... , y r +i) runs over X x V r with the condition 
rank(j/ 2 , ■ • • , Vr+i) + codim V = rank(x 2 , . . . , x r ) + codim 14 7 + 1. 

The theorem follows from this, because for 

fj,(r)= lim -log sup mult(/ t "|y)j /1 ,...^ r 

n^oo n . 

yieX,y 2 ...,y r eVcR d 

rank(y2,...,yr)+codim V—r—1 

we have: H sing (f t ) = fj,(l) < fi(2) < ■ ■ ■ < /j,(d + 1) = 0. 

To prove the proposition note that when r = we have (in this case we do 
not need W, V) : 

mult(/^) < (sup mult (ft)x) = f max sup mult (fj)x) ■ 

V£l ' \0<j<N xeX ' 

Let r > 1. Consider the continuity partition V Xl Xr for (ft) Xl ,...,x r , which 
is just the continuity partition Vx°}...,x r of {ft a )x 1 ,....x r (one iteration), and let 
^xt—xr f° r (/")xi,...,x r be the iterated partition. Note that the latter is the 
collection of all non-empty intersections 

P t0 n (f to r P l a (P* 1 ) n (f t0 )-} (/ tl )- t \ (P* 2 ) n • • • n (/^ . . . ) pt 1 „_ 1 (p 4 -), 

where P** are elements of Vxil..,x r and /p denotes the restriction of the (dif- 
ferential of the) map to the corresponding continuity domain. Every element 
of these partitions is invariant under the shift by vectors from span(x 2 , . . . , x r ). 
Therefore it intersects the unit sphere Si(x 2 , ■ ■ ■ ,x r ) in the orthogonal comple- 
ment. Consider the induced partition on the sphere and refine it so that every 
element has diameter no greater than e. Denote by n(V x ^ Xr , e) the minimal 

cardinality of such a refinement. Let also m(Vxil..,x r i £ ) be the maximal number 
of elements of V X \,...,x r that an e-ball B(y, e) D of S\(x2, • ■ • , Sr) 1 can meet. 

Denote by n(V fl W, e), m(V n W, e) the corresponding quantities in the 
subspace W. Then from the above formula for the iterated partition: 

n{V^x r nW,e)< n{V^...,x r n W, e) ■ m{V<fc+± n V, e), 

where Vl = f?( Xl ), y 2 = Vr = ft'M and V = ft'(W) with 

ft' = {R)xu..,x r - Therefore 

niK+'.'X n ^, e) < n(7>£* ... Xp n W, e) ■ sup sup m{V$£± n V, e), 

V yi,-..,y r 
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where the supremum is taken over all V C M. d and (t/i, . . . , y r ) € X X V Ar_1 such 
that codimension of (t/2 • • • , Dr) in V equals codimension of (x2 ■ ■ ■ , x r ) in 
Since for a fixed e the number n(V x ° e) is finite and 

yeS!(y u ...,y T )^nV 

the claim follows from the following statement. Fix i £ [0, N). 

Lemma 28. There exists e > (depending only on i) such that for all V C 
R d and all (yi, . . . , y r ) € X x V r ~ 1 , with rank(?/2, • ■ ■ , J/r) < dimV, and y £ 
Si(y2, ■ ■ ■ ,y r ) ± <^V there exists y' € (2/2, • • • , Vr) H V satisfying: 

\P}?L,Vr B (y> £ ) nV \^ ™*t((f i \v) Vl ,.,Vr,V')- 

This statement, modulo our notations and restrictions to V, is proved in |B2| . 
The proposition and hence the theorem follow. □ 

A. 3. Estimates on entropy by angular expansion rates 

It is possible to estimate the effect of angular expansion on the topological 
entropy of a piecewise afhne map / : X — -> X by its spherizations. Define the 
piecewise smooth map d% f ■ ST X X — * STft x \X given at x £ X' by the formula 



\dxf{v)\\ 



For x <E Sing(/) and v £ T x Sing(/) (the tangent cone) we let d x f(v) = 
lim.^d x s l ev f(v). For other (x,v) S STX the map is not defined. The angu- 
lar expansion of / is exactly the expansion in the fibers of its spherization. 

If d x f is degenerate we restrict to the orthogonal component of its kernel, 
and consider the map 

^/lKcr(d x /)-L : Ker(d I /) ± -> lm(d x f) . 

Then the map S x (f) = d x s) /| Kc r(d x /)-L : 5Ker(d x /) ± -> Slm(d x f) between 
(rank(d x /) — l)-dimensional spheres is given by the formula 

dx/kerfe/)-^) 



\\dxf\K C r(d x f)^(v)W 



For i < d we define 



pi(f) = lim - sup max log ||A fe d„S , :E (/ ,l )||. 

n^oo n ( XtV )0<k<l 

Let m* = muxr dimKer(d :E /) and d* = d — to* = max^ rank(d x /), where d is 
the dimension of X . The numbers Pi(f) can be non-zero only for i < d*. 
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We have: po(f) = 0. The number pi(f) measures the maximal exponential 
rate with which angles can increase under the map /. The numbers Pi(f) for 
i < d measure the maximal rate of expansion of the restrictions to i-dimensional 
spheres. If / is conformal, then p%{f) — for all i. 

Theorem 29 ([KR2 ). For piece-wise affine maps ff mu it(/) < J2t=i X Pi(f)- 
We define the maximal expansion rate 

AmaxCf) = lim sup -log || d x f n || , 

n—>oo x n 

and the minimal finite expansion rate 

Amin(/) = - lim SUp-loglK^/'^Kcr^/")^) -1 !! • 

In |KR2| we show that 

Piif) < » fA ma x (/) - A min (/)) . 

This gives the following result (the same bound holds for ht op (f)): 
Theorem 30. For a piecewise affine map f it holds: 

H sing (f) < A+(/) + d * (d * 2 ~ 1 \ \nUf) A min (/)) . 

B. Generalization of the Moran formula 

Consider an IFS (M d , /i, . . . , /jv) on a Riemannian manifold M, where the maps 
fi can possess singularities, but we assume that they are mild in a sense that 
the number of continuity domains is finite, any of them has piece-wise smooth 
boundary and the map, restricted to any of the domains, smoothly extends to 
the adjacent singularities (this is the case of the Zhang model). 

Remark that the IFS can be interpreted as the dynamical system (E^ x 
M,F), f(t,x) = (er^-t, ft x). Attractor y of the IFS can be defined via the 
attractor of the extended system F, which has the form A = E^ x y. 

Let || • || be the norm on TM generated by the metric on M . Denote 

s i = max 11^/,; ||, 8i = (maxH^/r 1 !!) . 

We assume that the maps are non-degenerate (this is just for simplicity of 
arguments) and strictly contracting, so that < < sf < 1. 

Let rj = max^={i | x — fi(jji) for some yi G y} be the maximal multiplicity 

of overlaps and = max#{y £ y \ x = fi(y)} be the maximal multiplicity of 

self-overlaps (we assume it is finite) on the attractor. Denote also by fii the 
multiplicity of the continuity partition for fi\y, i.e. the maximal number of 
continuity domains intersecting the attractor and meeting at one point of it. 
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Theorem 31. Let D_ = a, D = (3 be the solutions of the equations 

N N 

i=i i=i 

Then the Hausdorff dimension of the attractor satisfies: 

D < dim H (y) < D. 

In addition to Hausdorff dimension we will need some other dimensional 
characteristics (see |P2| for details). Denote by 9t(X, 8) the minimal cardinality 
of covers of X by balls of radius <5. Then the lower and upper box dimensions 
are defined by the formula: 

dim p (X) = lim log ^(X, <5)/ log i dlm s (X) = Hm log9I(X, 8)/ log ±. 

5^+0 ' .5^+0 

When these quantities are equal, their value is also called fractal dimension. 

Consider a Borel probability measure fi £ A4(X) (an SRB-measure on the 
attractor can be taken in the SOC-context if it exists). The upper and lower 
pointwisc dimensions are defined then as 

dJx) = lim log fi(B (x, 8))/ log 8, d ll (x) = lim log/j,(B(x,8))/logS. 

6^+0 8^+0 

When they are equal and constant a.e. the measure /x is called exact-dimensional. 
This is precisely the case, when supp fi = X and we have equality in the general 
chain of inequalities (together with (|13|1 below) : 

ess.inf d^x) < dim H pO < dim B (X) < dhn B (X), (12) 

where by essential infimum we mean its upper bound taken over all subsets 
U C X of measure 1 (and similar for ess. sup d fl (x)). The last two inequalities 
are known and the first one follows from the inequality ess. inf d „(x) < dimn(^) 
([EH), where dim H (^) = lim^+o inf{dim H (^) | fi(Z) > 1 - 8}. 

Other dimensional characteristics of the measure are defined similarly and 
satisfy: 

dhn B (u) < diras(M) — ess. sup ^(x). (13) 

Note that dims (/i) < dims (X), while the quantities ess. sup d M (x) anddims(X) 
are in general incomparable. 

The known formulas for the Hausdorff and other dimensions are generaliza- 
tions of Moran's result (|MIIH]) an d are based on the Bowen's equation (using 
the idea of coding) ; in this case one usually obtains exact-dimensionality |P2| . 
In the SOC-context coding becomes problematic in the presence of singularities 
(unless the properties of the SRB-measure are clarified) and thus we cannot 
easily establish exact-dimensionality or formula for the dimension. 

We prove instead the inequality of the theorem for all the various dimensions 
from l|12l) and (|13f) . which we denote just by dim(J^): 

D < dim(y) < D. 
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Proof. Let us consider at first the upper box dimension dims (J 7 ). The function 
91 satisfies the inequalities: 

The inequality from above is obtained as follows. Let S — {xj} be a <5-spanning 
set, i.e. a collection of points from X with Us(S) = X. Then fi(S) = {fi{xj)} 
may fail to be a 5 ■ s+-spanning set thanks to singularities. Whenever S <C 1, 
every (5-ball intersects maximally z9j domains of continuity for f t meeting y. 
Then we need to add maximally §i points for each ball U$(xj) intersecting 
singularities. The inequality from below is proved similarly. 
Now we have: 

y = fi(y)u---u.f N (y) 

and the same for [/^-neighbor hoods. This implies: 

N N 

9W) < E^W' 5 ) <J2^ y ,S/s+). (14) 



Denote a (6) — Vl(y , 5)5 dlmB ^' . This functions grows sub-polynomially: 

hm JSffW = 0. (15) 
<5^+0 log 1/5 

Lemma 32. Let Xi > 1 be some numbers and pi > be some probabilities, 
Yli=iPi — 1- Then \15\) implies: 



1 JV 



hm ^l^M < L 



Proof. Suppose the lower limit is > n > 1. Then for every sufficiently small 6 
there exists i € [1, N] such that a(XiS) > K,a(S). 

Denote A = max\<i<N Aj. Let C — max^-^ ^ a (5). Then: 

a{$) < -aiXiJ) < \a{X %1 XiJ) < ■< -4^-o"(A, 1 ...X ls(S) S), 

where s(<5) is the first number such that . . . Xi B(S) 6 G [1/A, 1]. This number 
can be estimated as follows: s(S) > — log 6/ log A — 1, whence: 

log a (8) logC — s(<5) logK log(CK) log ft 
log 1/5 ~ log 1/8 " log 1/5 ~ togT 

Therefore lima^+o 'log^/i — — isff < ^ an< ^ we S et a contradiction. This proves 
the lemma. □ 
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Now to obtain the inequality from above for dims (y) divide 114fl by yi(y, 5). 
Denoting zu = J2Zi^\ s t\ S ^ B{y \ *i = 1 /4 and P» = &i\s+\*^ B W /<& we get: 



' tr(XiS) 



1 X - 

w~^ P% o{8) 

1=1 



Thus Lemma |22 implies that 1/vj < 1 or 



N 



i<]>>i^i di 



imy 



and the first claim dims (y) < D follows from the contraction \sf\ < 1. The 
same arguments show another statement that c? M (x) < D. 
The inequality from below follows from 

N AT 

which implies $^ i=1 ^|Si~|^ lslB ^-' < ^ and the same for the lower pointwise 
dimension: d^(x) > -D a.e. 

In this case we should define 

(7(5) = yi{y,5)5^ {y) or o-(<5) =(ess.inf-log M (B(x, ( 5))) ( 5^^ 

respectively and use 

Lemma 33. Let Xi > 1 be some numbers and pi > be some probabilities, 
Sill Pi = 1- Then: 

lijs logaW =o ^ — EL^M> L 
log 1/6 s^+o a{5) 

This is proved similarly to Lemma 1321 The inequalities for the Hausdorff di- 
mension follows now from (|I2|I . □ 
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